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Abstract 


Doppler shift measurements of the Na D,; absorption line have revealed solar oscilla- 
tions in a new regime of frequency and wavenumber. Oscillations of vertical velocities 
in the temperature minimum and low chromosphere of the Sun are observed with fre- 
quencies ranging up to 9.5 mHz. There is no evidence for chromospheric modes of 
3 minute period. This indicates that the chromosphere does not form a good cavity 
for acoustic waves. 

The fundamental-modes appear with wavenumbers up to 5.57 Mm? (equivalent 
spherical harmonic degree, 3877). ‘The frequencies lhe below the predicted values at 
wavenumbers above | Mm7'. The values are in agreement with previous measure- 
ments that exist for wavenumbers up to 2.67 Mm7?. 

Spatial maps of velocity power show that high wavenumber oscillations are sup- 
pressed in active regions. ‘The shape of the power depression indicates that wave 
motion is affected in the laver of atmosphere where the measurement is made. The 
f-modes are suppressed in the same way as p-modes, indicating that the mechanism 
for wave suppression affects velocity fluctuations. Mode frequencies are not affected 


by the magnetic fields by more than 50 ;.Hz, the precision of the measurement. 
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Chapter 1 
Introduction 


Solar oscillations were first observed by Leighton (1960), and were explained by Ulrich 
(1970) and Leibacher & Stein (1971) as being the superposition of many global stand- 
ing waves. They have subsequently become an important tool for the determination 
solar structure. Measurements from helioseismology include the depth of the con- 
vection zone (Gough 1976), helium abundance in the convection zone (Christensen- 
Dalsgaard & Gough 1980), and internal rotation rate (Duvall et al. 1984; and others). 

In the past, observational limitations have restricted helioseismology to long wave- 
length oscillations (2 4 Mm). New instruments and techniques are beginning to make 
shorter wavelengths accessible to study. This opens the door to the exploration of 
structure at much smaller spatial scales. Stratification of the atmosphere, and vertical 
structure of magnetic regions are becoming active areas of investigation. 

The term solar oscillations conventionally refers to two types of waves. Because 
they exist in a resonant cavity, these waves are also called modes (eigenmodes of 
the cavity). Compressional or acoustic waves (p-modes) are longitudinal waves with 
pressure as the restoring force. They have at least one nodal surface in the radial 
direction. Gravity waves have gravity as the restoring force. A special case of these 
are the surface waves. Because they have no radial nodes, surface waves are called 
fundamental modes (f-modes). ‘These transverse waves propagate at the interface 
between two media. The remainder of the gravity waves are called g-modes. To date, 


g-modes have not been unambiguously detected, and they will not be considered here. 
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1.1 Wave Properties 


1.1.1 Acoustic Waves 


Sound waves can propagate vertically at frequencies greater than the acoustic cutoff 


or critical cutoff frequency w.. At lower frequency, a disturbance moves the whole 


atmosphere up and clown in phase. Hence, waves do not propagate. At low horizontal 


wavenumber & (A < si). the acoustic cutoff can be expressed in terms of the adia- — 


batic sound speed c¢ and the density scale height H. For an isothermal atmosphere, . 


the scale height is given by 


G : 
H =— (1.138 
Tg 
and, the acoustic cutoff, by : 
2 (1.2) 


where g is the acceleration due to gravity; and, y is the adiabatic index (Lamb 
1932). Taking c = 7.1 km/s. g = 278 m/s’, and y = 2 in the photosphere, we find 
H = 110 km and w,/27 = 5.2 mllz for & < 4.5 Mm”. 

In general, the acoustic cutoff depends on wavenumber. The relationship is given 


by (Kumar & Goldreich 1989) 


aan ae | Rae ee Jb 
ie (=—) 5 fll + akin le: 1+ kine — Beene (1.3) 


Figure 1.1 shows w, as a function of hk. 

Oscillations with wavelength greater than a few Mm propagate around the sun 
without excessive dissipation. A periodic boundary condition is, therefore, imposed 
on these waves. The wave equation is separable into radial and angular parts and has 


eigensolutions of the form, 
Waim(r.0.0.t) = Rulr) ¥"(8, 6) e (1.4) 


n=1,2,3.... (=0,1,2...79 SS eer 


Short wavelength oscillations are dissipated before traveling around the sun. More- 


over, at their small spatial scales, the medium can be approximated as plane-parallel. 





| 
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igure 1.1: Acoustic cutoff w,. vs. k 


Plane waves, therefore, provide a more natural basis set for the horizontal eigenfunc- 


tions. 
ey.) = Rax(r) eR Ft) (leo) 
The horizontal wavenumber / = |e] can be related to the spherical harmonic degree 
l as follows. 
(f+ 1) 
k= —————_ 1.6 
one ( 


Ray paths for a p-mode are diagrammed in figure 1.2. As waves travel downward, 
they encounter higher temperatures and an increasing sound speed. This causes them 
to refract back upward where they encounter increasing values of w,.. When the wave 
can no longer propagate (w < w.), it is reflected back down. The resonance condition 
is met when crests of downgoing waves are superposed on crests of refracted waves 
returning to the surface. For waves trapped in an adiabatic and polytropic layer 
(polytropic index = 4), this condition gives a dispersion relation of approximate 


form, 
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a~ 


é 


- 


Figure 1.2: Propagation of an acoustic wave in a plane-parallel atmosphere. The wave 
is refracted as it propagates downward, and reflected at the surface. The maximum 
depth reached by the wave is called the lower turning point. From Gough (1983). 


2 = oelat eons 
zs = (5+ g hye ee (a 
where g is the acceleration due to gravity at the surface. 


This dispersion relation can be measured directly from the ridges of power in a 
spectral decomposition (a plot of power vs. w and k) of solar velocities or intensities. 
The exact relation depends on the variation of sound speed with depth. 

The lower turning point d of the waves is approximately given by, 
3 (7 + 2) 
4k 


1 


(1.8) 


The excitation source of the p-modes is thought to be turbulent convection at 
the top of the convection zone. The efficiency of acoustic emission from convection 
is proportional to approximately the eighth power of the Mach number (Goldreich 
& Kumar 1988). Hence, only a few of the granules on the high velocity tail of the 
distribution are expected to contribute to the acoustic noise driving the p-modes 


(Brown 1991). 


1.1.2 Surface Waves 


Surface waves can propagate on the interface between convectively stable media (e.g., 
air over water). On the Sun, the transition region between chromosphere and corona 


provides such a boundary (Jones 1969). The model of Avrett (1992) describes the 


> 
— ———— 
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Figure 1.3: Temperature and density of the solar atmosphere. From Avrett (1992). 


temperature and density discontinuities shown in figure 1.3. For purely vertical den- 
sity stratification of a compressible fluid, surface waves obey the dispersion relation 
(Whitney 195s), 

Says (1.9) 


where g is the local acceleration due to gravity. 


The simple form of this relation has prompted the remark that, “There is no 
interesting physics in the /-modes.” It wasn’t until f-modes were found not to obey the 
above relation at high wavenumber that this conclusion was reexamined. Libbrecht, 
Woodard & Kaufman (1990) found f-mode frequencies to lie below the predicted 
values for wavenumbers above 1.4 Mm7'. At least three possible explanations exist, 
though numerical predictions have not yet been offered. 

Campbell & Roberts (1989) and Evans & Roberts (1990, 1991) predict that a 
horizontal magnetic field in the chromosphere would increase the frequencies of the 
f-modes by increasing the restoring force on the surface waves. The sign of the 


frequency shift may reverse at high wavenumber, however, in a more realistic model 
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of the chromosphere with an appropriate magnetic field structure (Roberts 1991). 
Rosenthal (1992) finds that evanescent waves dissipated in the corona would cause © 
power to leak out of the /-modes. This would lower the frequency of the modes at a 
given wavenumber by an amount dependent on the temperature difference between 
chromosphere and corona. 
Finally, structures in the atmosphere, such as flux tubes and convection cells, 
violate the assumption of horizontal homogeneity. It is not known whether this should 


increase or decrease the mode frequencies (Gough 1992). 


1.2 High Frequency Modes 


For the purpose of discussion here, the high frequency regime includes modes with 
w > w., where w,./27 = 5.2 mllz. This value of w, is the photospheric value for low 
wavenumbers. Because waves with w > uw, can propagate, they are not trapped by 
the photosphere. So, in the model illustrated by figure 1.2, there should be no upper 
reflecting surface at frequencies > 5.2 mHz. Removal of the top surface of the cavity 
eliminates the possibility of resonance, and hence, the mode structure. 

It was a surprise, therefore, when the p-mode ridges were found to continue above 
the acoustic cutoff (Jeffries ef al. 1988; Libbrecht 1988; Duvall et al. 1991). Two expla- 
nations have been put forward for the origin of the high frequency ridges. Balmforth 
& Gough (1990) propose that strong reflection from the base of the corona defines a 
chromospheric cavity which couples to the subphotospheric cavity. Waves with fre- 
quency high enough to pass through the photosphere would still be trapped by the 
corona. 

The presence of a chromospheric cavity has further consequences, however. The 
chromosphere is of order one wavelength in thickness for 5-minute modes; and waves 
of appropriate wavelength A, will resonate in this cavity. The p-modes propagate 
almost vertically in the atmosphere (angle of incidence ~ 0). Thus, 4, = c/y,, 
independent of horizontal wavenumber. The chromospheric modes should, therefore, 
have a dispersion relation. 


w = ears! (1.10) 
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Figure 1.4: Mode frequencies vs. wavenumber with a chromospheric cavity present. 


From Ulrich & Rhodes (1977). 


Ulrich & Rhodes (1977) found that the p-mode ridges would merge with a set of 
chromospheric modes at about 5.5 mHz. See figure 1.4. The p-modes would continue 
on above the 5.5 mHz band with an “avoided crossing” between the ridges above and 
below the chromospheric resonance. 

Dame, Gouttebroze & Malherbe (1984) claim to have seen the chromospheric 
modes in Ca II I\-line observations. Their power spectra are very noisy, however, and 
it is difficult to discern any definite ridges of power. Deming et al. (1986) find a peak 
in the one-dimensional power spectrum (power versus w) at 4.3 mHz. Its significance 
is difficult to assess, though, without knowledge of the wavenumber dependence of the 
power. Woodard & Libbrecht (1991) do not find this peak in similar one-dimensional 
power spectra. 

The other explanation for the high frequency ridges was offered by Kumar e¢ al. 
(1990). They propose that waves are emitted isotropically from a source near the top 
of the convection zone. Waves traveling directly toward an overlying observed surface 


interfere (in Fourier space) with waves refracted upward from below to produce “mock 
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modes.” It must be stressed that this interference is not physical. In fact, the directly 
emerging wave trains and refracted wave trains may arrive at the surface at different 


times and in different places. The basis functions, used to represent the waves, can 


only distinguish time, modulo 27/w, and position, modulo 27/k. Frequencies of the — 


mock modes can be used to determine the depth of the acoustic source layer (Kumar 
& Lu 1991). 
The two models for high frequency waves can be distinguished by observations 


that resolve the ridge structure in the 5.5 mHz frequency range. 


1.3. High Wavenumber Modes 


Most of the literature in helioseismology refers to the (spherical harmonic) degree 
! of a mode rather than the wavenumber &. By convention, low degree refers to 
0 <i < 4, intermediate degree, to 4 < / < 100, and high degree to 100 < l. The 
oscillations presented in this dissertation have 500 < 1 < 3877. By analogy with 
the electromagnetic spectrum, they should, perhaps, be called very high or ultrahigh 
degree, but I will stick to the term, high wavenumber. 

Previous observations of the /-modes exist up to wavenumbers of 2.67 Mm7! 
(Libbrecht, Woodard & Waufman 1990). As mentioned above, the frequencies do 


not agree with equation 1.9 at wavenumbers above 1.4 Mm7?. 


The variation of g 
with height accounts for only a 0.5% modification to equation 1.9. This is an order 
of magnitude less than the observed effect. [t would be desirable to confirm the 
observational result and determine the trend at higher wavenumbers. 

The size of structures that can be resolved by oscillations is roughly equal to 
the wavelength of the modes. For active regions, spatial scales range from 200 km 
(flux tubes) to 25 Mm (large sunspots). The corresponding wavenumbers range from 
0.25 to 30 Mm7’. Although the high end of this range is not attainable with any 
current instrument, there is still much room for investigation. 

Woods & Cram (1981) have reported a suppression of oscillation amplitudes in 
plage as compared to quict regions. Tarbell et al. (1988) and Title eé al. (1992) find a 


suppression of both convective motions and wave amplitudes. The suppression of the 






| 
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oscillations coincides with the boundary of the photospheric magnetic field to within 
one wavelength. 

Braun, Duvall & LaBonte (1987, 1988) and Braun et al. (1992b) have reported 
absorption of p-mode power by sunspots. The mechanism for this is not understood. 
Resonant absorption (Hollweg 1988; Chitre & Davila 1991), scattering to short wave- 
length (Bogdan & Zweibel 1987), and enhanced viscous dissipation (LaBonte & Ryu- 
tova 1992) are among the explanations offered. It is by no means clear whether 
only propagating acoustic waves below the photosphere are absorbed, or whether 
evanescent waves are also affected. A spatial map of oscillation amplitudes above 
the photosphere, where tle magnetic fields are different, would help untangle the two 
effects. 

Global oscillations have been found to have frequency shifts correlated with solar 
activity (Palle, Regulo & Roca Cortes 1989; Woodard et al. 1991). This may be due 
to a change in the shape of the cavity surface with sunspots present. Waves of a 
given frequency may have different wavenumbers within magnetic and nonmagnetic 
regions. Global measurements would then find a weighted average of the two disper- 
sion relations. Comparison of high resolution observations in quiet sun and active 
regions could resolve this issue by direct measurement of dispersion relations in the 


two types of regions. 


1.4 Height of Observations 


Oscillations can be observed as fluctuations in intensity and line-of-sight velocity. 
Both p-modes and f-modes produce motions of gas particles, and cause Doppler 
shifts of absorption lines. Pressure fluctuations due to p-modes produce tempera- 
ture/intensity fluctuations. The f-modes also produce intensity fluctuations, though 
the reason is not well understood. It may be that the ~ 100 km vertical displacements 
of parcels of gas affect the process of absorption line formation. 

Different parts of an absorption line are formed in different layers of the atmo- 
sphere. Hence, the height at which the velocities and intensities are determined 


depends on the positions in the absorption line at which measurements are made. In 
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Figure 1.5: Na D line profiles on left hand scale. Contribution functions for emergent 
intensity on right hand scale. Thick lines show the maximum of the contribution 
function at each wavelength. Thin lines show 1/3 of the maximum of the contribution 
function. From Schleicher (1976). 


the core of deep lines such as Na D,, very little continuum light emerges from the 
photosphere. The height of a measurement, therefore, depends in the contribution 
function for the emergent line emission (Gurtovenko, Sheminova & Sarychev 1991; 
Schleicher 1992). This function is shown in figure 1.5. 

The original intention of the observations described in chapter 4 was to make one 
set of measurements in tlie core and one in the wings of the Na D, line. Problems 
described in section 4.1 resulted in a measurement weighted mostly in the core (600- 
850 km), but with some contribution from below the temperature minimum (150- 


500 km). 
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Chapter 2 
Measurement Techniques 


Producing power spectra from uncalibrated images involves several stages of data 
analysis and reduction. ‘|This chapter reviews some general procedures, standard and 
nonstandard, for measuring line-of-sight components of velocity and magnetic field, 
and for determining wave amplitudes and phases. The details of calibration to remove 
artifacts of the instrument and atmosphere are left to chapter 5. 

The line-of-sight velocity is found from the Doppler shift of an absorption line. 
The Doppler shift, in turn, is determined by measuring intensity as a function of 
wavelength at a series of fixed-wavelength positions in the absorption line. The more 
points in the measurement, the better the velocity estimate will be. Unfortunately, 
more points require proportionately more time to make the measurements, and geo- 
metrically more time to analyze them. 

In this chapter are described two new techniques for making velocity images (dopp- 
lergrams) and line-center intensity measurements using three-point and four-point 
methods. The four-point method was developed for the Solar Oscillations Investi- 
gation / Michelson Doppler Imager (SOL/MDI) experiment for the Solar and Helvo- 
spheric Observatory (SOHO) spacecraft. It was not applied to the data used for this 
dissertation. 

Images of line-of-sight magnetic field strength (magnetograms) are an important 
tool in this study. The tried and tested “two-point magnetogram” method is described 


below. 
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Figure 2.1: Measurements of the Na D, line 


Finally, the Fourier transform is essential to the study of waves propagating in 
a linear medium. Procedures used to extract information about the oscillations are 


described in the last section of this chapter. 


2.1 ‘Three-Point Measurements of Line Parame- 


ters 


The three-point method for finding Doppler shifts and line center intensities assumes 
a parabolic absorption line profile. This assumption is valid when the three intensity 
measurements are positioned relatively close to the line core. Figure 2.1 shows the 
measurements schematically for the Na D, (5896 A) line. The Doppler shift and line 


center intensity are the r and y coordinates of the vertex, respectively. 
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2.1.1 Algebraic Solution 


Assuming a parabolic line profile 


HO = ee (A=) (2.1) 


const, 


we can solve for the vertex (Ao. /p) in terms of three intensities /,, /2, /,; measured at 


wavelengths A;, A. A3. respectively. 


os BA + A3C3 


Cie ca 


(dc, = 63¢3)° 


io =) 9 SS oe) 
ae BO We aie C10) (23 C1) 28) 
where a = (Ay + Ae 2). Ce = a(A3 =~ Ao), 
Oy = Ay = AQ. 63 = A3 — Ao. 


an = 0, ( 13 = I>), and c. = 63( 1; = Ip). 


2.1.2 Error Estimates 


The errors on the measured quantities /,./,./3; are due to limitations of the instru- 
ment. The transmission of the tunable filter varies with wavelength as shown in 
figure 3.8. The phase of each Fourier component of this modulation is independent, 
and drifts with time. In addition. the wavelengths tuned by the filter are quantized 
in units of the smallest step size (3.14 mA). These sources of error were used in a 
simulation to estimate errors for the velocity and intensity. 

The simulation assumed a parabolic line profile with parameters measured from 
the data. A range of velocities of +2000 m/s with a flat distribution, and a range of 
line-center intensities of £8% with a Gaussian distribution were used in the model. 
Intensity “measurements /,. /;, /; were modeled assuming a random jitter of up to 
+1.57 mA (flat distribution) in the location of each measurement, and assuming 
phases for the transmission modulation which varied randomly between sets of 3 
filtergrams, but not within a set. The velocities and intensities were then calculated 


using equations 2.2 and 2.3. ‘The errors had approximately Gaussian distributions 
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with standard deviations. 


= 300 m/s. anda — (sae (2.4) | 


2.2 Four-Point Measurements of Line Parame- 
ters 


An algorithm was developed for the SOI/MDI experiment to determine Doppler shifts 
using four measured intensities. The measurements are used to solve a system of si- 
multaneous equations for the four absorption line parameters. Since computing time — 
must be kept to an absolute minimum, the nonlinear equations are solved (numeri- 
cally) using a look-up table approach. 

The model included ouly the Lyot and Michelson filter elements, and the following 


parameters were used: 


e Lyot filter: profile - - Gaussian 


fwhm -- 500 mA (see section 2.2.3) 


e Michelson elements (2): profile — cos?(kA) cos*(2k\) 


[voip eels ae 
4 iL 


e Nil (6768 A) Line: profile — Gaussian 
fwhm 130 mA 
depth — 58% 


e Measurement Spacing — 85 mA (symmetric about nominal line center) 


See figure 2.2 for illustration. 
Four measurements are made at increasingly longer wavelengths across the ab- 
sorption line. to yield intensities /; through /,. The ratios, a and 3, are formed as 


shown: 


ere 
Fisk 
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Figure 2.2: Spectral profiles. 
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The quantities a and 3 are independent of CCD dark current and gain variations. 


(2.6) 


They are also independent of continuum intensity and line-depth if a suitable choice 
is made for the Lyot filter profile (see section 2.2.3). 

We can solve exactly for the Doppler shift Ao, line depth d, and continuum inten- 
sity [9 if we assume a cosine shaped absorption line with measurement spacing equal 


to the hwhm: 


ON) = he [ - : (co SE — i) (2.7) 

2 whm 

Solving, we find, 
Ao = fwhm - (; = Goi! a (2.8) 
l 

cs all +], —s) (2.9) 
-— 2.10 
i= 2 (2.10) 


where s = ,/(/, — [3)? + (14 — I)? 
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Figure 2.3: Ratio a vs. velocity for width increments: 0, 10, 20, 30%. 


mm 


Note that @ has a one-to-one correspondence with the ratio r = () 5) =e 
(41 —13 )+ (4g = 12) 


which appears in the SOI/MDI proposal (i.e., r = (a—1)/(a+1)). Hence determining 
the velocity from a alone is equivalent to determining it from tan7!r. 


For a Gaussian line-profile, 
Lalas 
hOotauk f eee (2.11) 


Where o = fwhm/2.36 


figures 2.3 and 2.4, respectively. It is apparent that the velocity determined from 
a alone would be in error if the calibration table were not made for the “right” 
line-width. 

This caveat can be removed by determining the velocity v from both a@ and 2: 


| 
| 
The ratios a and @ are plotted as functions of velocity for different line-widths in 
v = v(a, 3). An algorithm for doing this is outlined below. 


CHAPTER 2. MEASUREMENT TECHNIQUES 17 
2.0 


0.8 





—1000 0 1000 2000 £3000 4000 5000 
velocity (m/s) 


Figure 2.4: Ratio 3 vs. velocity for width increments: 0, 10, 20, 30%. 


2.2.1 Implementing the Algorithm 
Look-Up Table Approach 


It is possible in principle to compute velocity by using a two-dimensional look-up table 
in a and BG (ie., v = v(a,Z)). This, however, would require an inordinate amount 
of memory ((12 bits in a) x (9 bits in 8) = 2M words) to achieve an accuracy of a 
few m/s. A better strategy is to take advantage of the fact that a, alone, provides a 


pretty good estimate of velocity; and write: 
v(a, 3) = vy(a) + v(a, 8) (2.12) 


where v2 is a correction to the velocity v, and, hence, does not vary with a@ as quickly 
as does v}. 

This approach can achieve an accuracy similar to the previous approach but with 
much less memory overhead. The table, v,(@), requires a 12 bit address; and the 
table, v2(a, 2), requires 7 bits in a by 9 bits in @ (this has not been optimized); for 
a total of 66k words. 
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Domain and Range of a and 3 


The ratio a has a pole in the negative velocity-domain (i.e., a red-shifted line). To 

avoid this. a test is made to determine the sign of the velocity. For positive velocities: 
a = 4-4, For negative velocities: a = 4-8, 
=e hy 

The velocity is negative if and only if: 4, — 13 > I, — I, (2138 

As can be seen from figure 2.3, @ covers a more-or-less symmetric range about 


zero. Therefore. 1t can be scaled-up to provide an address for a look-up table without 





having to subtract an offset. 
The function 3 is symmetric in the velocity domain and is continuous for sym- 
metric line profiles. For asymmetric lines (and noisy lines) 3 has a pole near v = 0. 
In practice this does not cause problems because. near v = 0, v is independent of 8 
(1.e., vo(a, 3) & 0 for small vc). An offset may have to be subtracted from 8 when 
scaling it for use as an address. 
Finally, note that: 
elk & Seetias  Waamle elle ieee’, 


bps Se 
ieee (i a 


=| (2.14) 


Hence, both a and 3 can be obtained from the differences (/; — /3), (4; — Jy), and 


CMe) 


The Algorithm 


It is assumed here that tle differences: A, = 1; — I5, A> = J; — la, and Age eee 
have been accumulated in image buffers. 


For all pixels: 


Calculate 3 and sign of velocity: 


numerator = A, — +43 
| —l. if numerator > 0: 
sgn = 


+1. otherwise. 


numerator “4 
Pe 
A> l 
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Calculate a: 
oe it gna: 


aad as 
= es 
Ay 


<2. otherwise. 
Calculate velocity: 


v= sign -[V4(C2a) + 12(C3 a, Cy 2)} 


C,,...,Cys are constants: 1) and V2 are look-up tables. 
It should be noted that noise in the data may cause a “division by zero” error in 


the algorithm. 


Constructing the Calibration Tables 


The following procedures are used to construct the calibration tables , and V. 
Parameters used for this study appear in square brackets. 


The table (a) is made for a line-profile of nominal width {130 mA] as follows: 


1. Calculate a for a representative number [44] of discreet (evenly spaced) 
velocities ranging between some number slightly less than zero [-40 m/s] 
and some number shehtly larger than the largest velocity expected from 


the data [4500 m/s]. 
2. Invert the table of a versus v to get a table of v versus a. 


3. Interpolate (cubic spline) between the endpoints to get V,(q@), a table of 


velocities for an evenly spaced abscissa a. 


The table V3(a,) is made for the same range of velocities as is Vj; and, for a 
given range of widths [100 to 295 mA]. The minimum width should not be less the 
measurement spacing [85 mA]. 

There are two “unphysical” regions in the table 14. The highest values of a cannot 
occur with the highest values of 3, and the lowest values of a cannot occur with the 
lowest values of 3. These portions of the table will not normally be addressed by the 
Doppler algorithm. 


VY. is calculated as follows: 
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1. Calculate a and 3 for a representative grid of velocities [22] and widths 
[14]. 


2. Interpolate (spline) to get a much finer grid [512 in velocity by 256 in 
width]. 
3. For each “real” velocity v (512 of them] calculate the velocity correction 
os 
vg = vu —V{(a) 
using the table })(a) computed above. 


4. Sort the ordered triplets (v2, a,.3) into a 2-dimensional table Va(a, 3). keep- 


ing track of how many times each bin gets filled. 


ay | 
* 


Fix uncooperative bins: 


a) Ifa bin accumulated more than one value of vg, normalize to get 





the “average” vy for that (a, 2). 

b) Ifa bin in a “physical” region did not get filled [there were 6 of 
them]. interpolate [linearly] using adjacent bins. 

c) Fill sunphysical” bins with the value of v2 from the nearest “physi-— 
cal” bin. This is done because noise in a measurement may produce 
an ordered pair (a,.3) which does not le quite where it should in 


the talsle: 


Tables 1; and \3 were made in units of m/s, and rounded to the nearest integer. 


2.2-2 Spectral Simmiation 


The instrument and Doppler algorithm were simulated in order to determine the 
magnitude of statistical and systematic errors. The simulation program stored the 
line and transmission profiles in single-precision 512-element tables covering the range 
+800 mA about nominal line-center. The intensities. I,....,24, were calculated tones 


discrete set of velocities as follows: 


1. Offset the index of the line-profile table to simulate a Doppler shift. 
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2. Take the product of the transmission-profile and the offset line-profile. 


3. Integrate [sum] over wavelength to get “measured” intensity. 


The Riemann interval and limits of integration were varied to verify that they 


did not affect the results. 


The error estimates are described in the next two sections. 


2.2.3 Sources of Systematic Error 
Line-Depth Variations 


The ratio a will be sensitive to variations in line depth if measurements /; and J, (or 
I, and /,) do not integrate the same amount of continuum. 

Figure 2.5 shows the transmission through the Lyot and Michelson filters as a 
function of wavelength for the first and third measurements. The ratio of integrals of 
these curves 1s plotted as a function of Lyot filter fwhm (assuming a Gaussian Lyot- 
profile) in figure 2.6. The two integrals are equal when the Lyot fwhm is ~ 625 mA. 
They are different by 2% when the fwhm is 500 mA. 

This 2% difference causes a velocity error €,which depends on line-depth as shown 


in figure 2.7. If the Lyot fwhim were 400 mA, €, would be larger by a factor of 4. 


Lyot Filter Drift 


If the Lyot filter profile drifts in wavelength space, the four measured intensities may 
not integrate the same amount of continuum. This would give erroneous values for a 
and 2 and, hence, v. To eliminate this source of error the transmission profile ti, 5:() 


should satisfy: 


miich 


J tye ~ 6) eee | ered (2.15) 


for some small range 6, over which the Lyot filter might drift; t!,,.,(A) 1s the trans- 
mission of the Michelson filter at wavelength setting 7: 
For a Gaussian Lyot-profile. it turns out, this condition is more or less satisfied 


for the same fwhm that minimizes ¢, due to line-depth variations. The reason can be 
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Figure 2.5: Transmission product for Lyot and Michelson filters vs. wavelength. ¢, 
and t3 are the transmission profiles for measurements /, and Js. 
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Figure 2.6: Transmission ratio f t3dA/ ft; d\ as a function of fwhm of the Lyot filter. 
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Figure 2.7: Systematic error €, vs. velocity for line-depth increments 


—20, —10,0, 10, 20%. 


seen from figure 2.5. For transmission profile, t,, if the Lyot profile shifts to the left 
(right), any increase (decrease) in the area under the main-lobe is compensated by a 
decrease (increase) in the area under the largest side-lobe. Transmission profile t3 is 
unchanged to first order since it lies close to the maximum of the Lyot-profile. 


Figure 2.8 shows ¢, for 8 mA offsets of the Lyot filter. 


Line-Width Variations 


As stated above, if a alone is used to compute v, an error will be introduced which 
depends on the line-width. Figure 2.9 shows the error incurred if the line-width is 10, 
20, and 30% larger than the nominal width for which the calibration table was made. 


If both a and @ are used. the error is reduced to a few m/s as shown in figure 2.10. 
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Figure 2.8: Systematic error €, vs. velocity for Lyot drifts: —8,0,8 mA. 
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Figure 2.9: Systematic error €, vs. velocity for line-width increments: 0, 10, 20, 30% 
(a-only method). 
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Figure 2.10: Systematic error €, vs. velocity for line-width increments: 0, 10, 20, 30% 


(a, 8 method). 


Line-Asymmetries 


The effect of line-asymmetry was explored by using a line-profile with a Gaussian- 
shaped bisector as shown in figure 2.11. It was found that the depth at which the 
velocity is measured depends on the velocity, but is independent of bisector shape. 
The velocity dependence is plotted for different bisector shapes in figure 2.12 (a-only 
method), and figure 2.12 (a. 3 method). Note that ¢€, does depend on the precise 
shape of the bisector, but its magnitude cannot be determined until a measurement 


is made of the bisector shape at the spatial resolution of the MDI instrument. 


2.2.4 Statistical Error 


Statistical errors were found by a standard propagation of error calculation. Deriva- 


tives of the functions v;(a) and v2(a. 3) were approximated by finite differences. 
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Figure 2.11: Geometry of modeled line asymmetry. 
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Figure 2.12: Depth of measured velocity vs. velocity for bisector shape hwhm: 
=e 7 a (a-only method). 
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Figure 2.13: Depth of measured velocity vs. velocity for bisector shape hwhm: 
—7.5,...,7.5 mA (a. 3 method) 


For a sinusoidal line-profile with delta-function transmission-profiles, the statist:- 


cal error o, can be written down in a form like that shown in the SOI/MDI proposal. 


Wesel 2 — a 
iy [one 


where, s is the measurement spacing; d, the line depth; and, Jp, the number of photons 





Oo, = 


(2.16) 


collected at continuum intensity (half the well depth for the CCD). 

Numerically, ¢, from this calculation is 16 m/s (with d = .58, Ig = 10°). The finite 
spectral resolution, however, tends to wash out the signal. For Gaussian transmission- 
profiles (fwhm = 100 mA), o, is 25 m/s. For the transmission-profiles in figure 2.5, 
a, is 30 m/s as shown in figure 2.14 (a@-only method) and figure 2.15 (a, @ method). 
The scatter in figure 2.15 is due to the method by which the calibration table V2(a, 3) 
is calculated. 

The sudden increase of o, for velocities above 3500 m/s is due to the increasing 
dependence of v on J together with the fact that 8 is only allotted 9 bits in the 


calibration-table address. Allotting more bits to @ would proportionately reduce this 
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Figure 2.14: Statistical error a, vs. velocity for different line-widths (a-only method). 
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Figure 2.15: Statistical error a, vs. velocity for different line-widths (a, 8 method). 
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behavior. 


2.2.5 Discussion of the Method 


The studies indicate that velocity errors due to line-depth variations and Lyot filter 
drift can be minimized by proper design of the Lyot filter profile. In principle it 
should be possible to eliminate the errors entirely, though only a crude optimization 
has been done here. 

Line-width variations can produce velocity errors of up to a few hundred m/s 
when the velocity is calculated from a single rational function of the data (a-only 
method). The error is reduced to a few m/s by using two rational functions (a, 2 
method). The more accurate computation requires an extra image buffer and entails 
doing 2.2 times as many operations. 

A line-asymmetry will affect the measured velocity because the line-bisector is 
sampled at different depths across the solar disk. This happens because the average 
velocity is a function of disk position. and the data tend to sample more of the 
line-wings for large Doppler shifts. 

Data taken at La Palina in 1988 show that the averaging of subresolution velocity 
elements produces a “noise” level of ~ +20 m/s. Hence, it seems reasonable to set 
this as an upper limit for systematic errors. 

The statistical error for a four point measurement is ~ +30 m/s with a CCD 
well-depth of 2 x 10° photons. The error becomes significantly larger at v ® 3500 
m/s for the a, 3 method. This can be rectified by use of a larger calibration table. 
The errors for 2 pass (8 point) and + pass (16 point) measurements are 21 m/s and 
15 m/s, respectively. 

In order to minimize the statistical error, sidelobes should be kept to a minimum. 
Hence the optimum Lyot filter-profile would have a flat peak with wings just large 
enough to ensure integration of equal amounts of continuum at each wavelength 
setting. 

This recommendation for the Lyot filter width has been used in the design of the 
MDI instrument. The choice of velocity calculation method is still pending further 


study. It is not yet known whether lne-width or line-asymmetry variations will 
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contribute more to measurement errors. Only if line-asymmetry variations are small” 


can the a, 2 method used. 


2.3 Magnetograms 


2.3.1 Inverse Zeeman Effect 


Spectral lines of atoms in a magnetic field are split by the Zeeman effect. For a field | 
component fy pointing toward the observer, an emission line of wavelength Xo is split | 
into two components (go components) with opposite circular polarizations. The left | 
and right circular polarizations (LCP and RCP) are displaced from Ag by —A) and | 
+A\, respectively. 


ine | 
Ay = — Ee — 9 (2.17) 


“ Srm.c? 
where g’ and g” are the Lande factors. and VW‘ and MM" are the magnetic quantum | 


numbers for the initial and final states (Bray & Loughhead 1964). The magnetic } 


quantum numbers are subject to the selection rule Wf" — Af’ = +1. For the Fe @ 
(6302.508 A) line. 
eX . 
ra = 17.4 mA/kgauss. (2.18) 
‘Il 


For absorption lines. the LCP and RCP components will be preferentially absorbed 


at Ag — Ay and Ap + A), respectively. This is termed the inverse Zeeman effect. 


2.3.2 Longitudinal Field Measurement 


From figure 2.16 we can write down the measured LCP and RCP intensities. 








[ 
= 2.19) | 
pias ||| (2.19) | 
Ol 
As } ze 





where - is assumed to be constant over the range \,, + A). 





CHAPTER 2. MEASUREMENT TECHNIQUES 31 


om, 
& Cad 


ue RCP — 


» 


A =o uae ASX , 
AN, N 
: 


Figure 2.16: Measurement of Zeeman splitting 


We can now define the magnetogram signal .V, 


ie Tat Ol 








= = — |—/A Pa 
I, + Ir Io OX ; | ) 

Combining this with equation 2.18, we have, 
Ay =cM (2.22) 


where the constant of proportionality c depends on the profile of the absorption 
line and of the filter through which the intensities are measured. The constant is 


independent of continuum intensity. For the filter parameters in table 4.1, 


bo 
bo 
W 
See” 


eS 0 eauss: (2: 


2.4 Transform Methods 


The goal of the data analysis is to measure the dispersion relation and amplitude 


of waves. The natural method for plane wave decomposition in a linear medium is 
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the Fourier transform. Consider a monochromatic plane wave p with wave vector — 


ko = (Nore ha) aonlaltice cqueniia aresny 


a? 


p(v.t) = cos(ko - ¥ — wot) (2.24) 
The Fourier transform of p is defined as, 


Fe {p} = [ / a cos(kg » F — wot ei Ft 2x dt (2.25) 


Evaluating the integral. 


ae eee ere ae 
Fi =0°(K — hg) O(w — wo) + 50 (A + ko) 6(w + wo) (2.26) 


— om 


The two delta functions, lving opposite one-another in (k,w) space, illustrate a general 
property of Fourier transforms. For real valued functions, the (complex) Fourier | 
transform will have a factor of 2 redundancy in information. In fact, we only need 
to look at a half-space (containing the origin) to obtain all the information about 
the original function. It is conceptually easiest to eliminate the negative frequency 


half-space, and consider only: 
=O hy OC. =o <h,< 0, 0 Sage 


It is important to account for relative motion of the observers field of view with | 
respect to the observed medium. If the medium is moving at constant velocity v with 
respect to the observer, we have * — x — vt. This produces a measured frequency | 


shift w — w+ A, for the plane wave. where 
A) = aa (2.278 


A final and extensively used property of Fourier transforms is the convolution 
theorem. It states that the Fourier transform of the convolution of two functions — 


f(x) and g(x) is equal to the product of the Fourier transforms of the functions. 


lo 
CO 
~~ 


Fi f f(g) gle —Gde} = FUS(@)} Fala) (2. 


Xx 
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2.4.1 Practical Considerations 


The discrete Fourier transform. or Fourier series. can be viewed as an extreme example 
of a Fourier transform (Bracewell 1986). The series representation F, is useful for 


discrete, evenly spaced data f, with periodic boundaries. It 1s defined by, 


: I = —2imnk/N 
fe =Fi{fas= Tv ene ue e723) 


| 
The inverse transform is defined by. 


l ASS ae a 
iz == jae i = = ; Pee" nk/[N (2.30) 
oe) 


"I 
For measurements mace at an interval 6. the maximum Fourier frequency that can 


be represented is the Nyquist frequency. 
| l se 
lee a ay) (eon 


The discrete transform has the advantage that it can be evaluated using a very 
efficient algorithm. the Fast Fourier Transform (FFT) algorithm. The data from the 
series of Doppler and intensity unages is indeed binned discretely on an evenly spaced 
grid in the spatial dimensions (by pixel) and temporal dimension (by image number). 
The condition of periodic boundaries must be enforced, however, by apodizing (see 
below). 

There are several algorithms available for discrete Fourier and related transforms. 
For handling real valued data. it 1s more efficient (by a factor of 2) in computing 
time and memory usage to construct the Fourier transform from a Fast Sine-Cosine 
Transform. For doing convolutions to smooth and filter data, it is most efficient and 
conceptually simplest to use the Fast Hartley Transform. The abbreviation FFT is 
used generically hereafter to refer to the Sine-Cosine and Hartley transforms. 

The FFT algorithms used were designed for data series containing N points, where 
N is even and factorable into low primes. The spatial dimensions were 512 pixels in 
both z and y (.V, = .V, = 2”). The temporal dimension was 480 images (.V,; = 2°-3-5). 
This number was chosen because it entailed the largest “nice” number of filtergrams 


that would fit on a single data tape. 
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2.4.2 Apodization 


To approximate periodic boundaries, the data is multiplied in each dimension by a 
function which tapers to zero at the edges. The exact form of the function must 
be chosen very carefully to avoid disturbing the Fourier spectrum. The convolution 
theorem tells us that multiplving the data by an apodizing function is equivalent to 
convolving the Fourier spectrum with the Fourier transform of the apodizing function. 
Therefore, the apodizing function should be constructed such that its transform has 
a narrow bandwidth and very small sidelobes. 

Several standard apodizing functions are listed in Harris (1978). The 3-term — 
Blackman-Harris function was chosen. Its discrete Fourier transform has a 3 dB | 
width of 1.56 bins. and largest sidelobes of —61 dB. The windowing function for N 


bins is given by 


2s BU 
lV =" er, COS (Sn) + (2 COS (=) (2.32) 
A ay 
A ke a N 
{ol t= == 5—— 7 lose ee 


where dg = 0.44959 

a, = 0.49364 

i =O ors 
When filtering data as done in chapter 7, it is not, in principle, necessary to 
apodize the data. For the La Palma data sets. however, the tmage rotates during 
the observations. Removal of this rotation from the images leaves the corners of the 
images without real information. To minimize round off errors from the corners, the 
data were masked before filtering. The mask used had a cosine taper (to zero) at the 


edges of a flat central disk: 


| if ronn < iN ae ee 
oan i U if Pmn = Nv 
cos (Sinn — N+ 32)) otherwise 
wihtere 7,4 “= oleae 
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2.4.3 Normalization 


The FFT’s are normalized such that the spectral power of a sine wave 1s equal to the 


square of the wave-height before apodization. 


1 


FL, > sin kon \| == | (2.34) 
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The exact normalization is a weak function of & and is different for the real and 


imaginary parts of the transform. However, it can be approximated by a constant. 


IF {fa} F = AIF fa} | (2.35) 
. il 
Nay 
where A= V7.) (2.36) 
n=- = 


This approximation is valid for wavelengths that are short compared to the scale 
of variation of the apodizing function JV). as is the case for all wavelengths and 
frequencies of interest. 

There is an additional wavenumber-dependent normalization factor that corrects 
for the attenuation of spatial frequencies by the atmosphere and optics. This correc- 


ments cescribed in sections 3.1 and 5.4.1. 


2.4.4 Statistical Error Estimates 


An error estimate can be obtained for the spectral power |F,|? by a standard propa- 
gation of error calculation. If we assume a constant error of on the data f,, we can 


use equation 2.29 to find the error on the power oj, ;2. 
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Chapter 3 


Telescope and Instrumentation 


3.1 Telescope 


The 50 cm Vacuum Solar Telescope at the Swedish Solar Observatory is shown 
schematically in figure 3.1. A photograph of the tower at Roque De Los Mucha- 
chos on the island of San Miguel De La Palma (Canary Islands, Spain) appears in 
figure 3.2. The telescope is a refractor with an f number of 45. 

The image forming lens is a 50 cm achromatic doublet. Two flat mirrors M1 and 
M2 in the altazimuth turret guide the beam down an evacuated tube. The first lens 
of the achromat forms the vacuum window. At the bottom, the flat mirror M3 directs 
the beam through a vacuum window to the (Ursies) observing room. Extreme care 
was taken by the designers of the optics to ensure the telescope would perform at the 


diffraction limit (Scharmer et a/. 1985). 


3.2 Postfocus Instrumentation 


3.2.1 Tunable Filter 
Lyot Filter Operation 


A Lyot filter is constructed from a series of simple birefringent elements. A represen- 


tation of one such element is shown in figure 3.3. An incident pulse of light is polarized 
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Figure 3.1: Schematic diagram of the Swedish 50 cm Vacuum Solar Telescope. From 
Scharmer et al. (1985). 


at 45° to the optical axes of a birefringent crystal. Fast and slow pulses emerge from 
the crystal, and pass through an exit polarizer. The output of the element is, thus, 
a pair of identically polarized. equal amplitude pulses. with a temporal separation 
proportional to the crystal thickness. The spectral response of such an element is 
given by the Fourier transform of its impulse response. The Fourier transform of a 
pair of delta functions is a cosine. 

A Lyot filter consists of a cascade of elements, each double the length of the 
previous one (see figure 3.4). The transmission of each element as a function of 
frequency is shown in figure 3.5. The longest element has the most closely spaced 
maxima. The next longest has minima that fall on every second maximum of the 
longest element. and so forth. The net effect is to produce a transmission profile 
with a mainlobe width approximately equal to that for the longest element; and a 
free spectral range. or distance between peaks, equal to the period for the shortest 


element. 
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Figure 3.2: Swedish Solar Observatory, La Palma, Spain. 





CO 


CHAPTER 3. TELESCOPE AND INS PRU ion ties 39 


INPUT 
PULSE 
ENTRANCE POLARIZER 





Figure 3.3: Pulse propagation through a simple element. From Title & Rosenberg 
(1981) 





Figure 3.4: Lyot filter with four elements. From Title & Rosenberg (1981). 
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Figure 3.5: Transmission vs. frequency for the components of a Lyot filter (A, B, 
C, D) and the transmission through successive stages of the cascade (E, F, G). From 


Title & Rosenberg (1981). 


CHAPTER 3. TELESCOPE ANDUINSTRUVE ie 


EXIT POLARIZER (0°) 


HALF WAVE PLATE 


ON 


(ROTATES) 


QUARTER WAVE 
PLATE (0°) 


BIREFRINGENT PLATE (45°) 


ENTRANCE POLARIZER (0°) 





Figure 3.6: Birefringent tuning element. From Title & Rosenberg (1981). 


Tuning a Birefringent Element 


| 
. 


The bandpass of a birefringent filter can be tuned over the full free spectral range 


by individually tuning the component elements. The best practical method for doing 


this is illustrated in figure 3.6 (Title & Rosenberg 1981). As in figure 3.3, a fast 


and slow pulse emerge from the crystal parallel to the fast and slow axes. The delay 


between the pulses is inversely proportional to wavelength. Reconstituting the two 


polarizations gives elliptically polarized light. The axes of the ellipse are horizontal 


and vertical, and the ratio of the axes, the ellipticity, is a function of wavelength. 


A quarter wave plate transforms the polarization from elliptical to linear, and the 


orientation of the resulting polarization vector is dependent on ellipticity. Finally, 


a rotatable half wave plate and fixed exit polarizer are used to select the linear 


polarization, and hence, the wavelength. 
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The Solar Optical Universal Polarimeter 


The Solar Optical Universal Polarimeter (SOUP), developed by LPARL, employs a 
nine element tunable birefringent filter. It 1s stmilar in design to a Lyot filter, but each 
birefringent element contains a partial polarizer to apodize its transfer function. The 
result is a filter with smaller sidelobes than a Lvot filter. The SOUP instrument was 
used on the Spacelab 2 shuttle flight, and the filter has become the evaluation model 
for the Orbiting Solar Laboratory — Coordinated Instrument Package (OSL-CIP). 
Software-selectable blocking filters. each 6-8 A wide, are centered on solar absorption 
lines. Waveplates can also be selected to measure left and right circular polarizations 
as well as linear polarizations spaced every 45°. 

At a wavelength of 5896 .\. the free spectral range of the tunable filter is 14 A, 
and the fwhm of the passband is 62 mA. A “wide mode” is also available by rotating a 
half wave plate within the longest element. This effectively shortens the element by a 
factor of 3, and increases the bandwidth to 103 mA. The stepping motors that rotate 
the tuning half wave plates permit tuning in 3.14 mA steps. Tuning is repeatable to 
within half a step. 

Temperature gradients and imperfections in the filter components produce small, 
but measurable, variations in filter characteristics over the field of view. The center 
of the passband varies by up to 10 mA and the fwhm, by up to 5 mA over the field of 
view. In addition, the total amount of light transmitted (integrated over wavelength) 
is modulated as the filter is tuned. (See figure 3.8.) This modulation accounts for 


most of the error in the velocity and intensity measurements. 


mee.2 CCD Camera 


The CCD (charge coupled device) camera uses a Texas Instruments 1024 x 1024 pixel, 
virtual phase detector. It is also an evaluation model for OSL-CIP. Subareas of the 
full detector can be read out to reduce the time required to write data to tape. In 
addition, 2 x 2 groups of pixels can be summed to decrease exposure time by a factor 
of 4. This 2 x 2 summing mode was used for the majority of observations listed in 


chapter 4. 
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Figure 3.7: Transmission profile of the SOUP tunable filter vs. wavelength. Trans- 
mission was measured using a spectrometer, with the filter tuned to 6328 A. 


The detector geometry consists of a rectangular grid of pixels. Each pixel is 
efficient at collecting light up to the boundary with its adjacent pixels. Thus, the full 
area of the CCD array is sensitive to light. 

The charge accumulated on a pixel is proportional to the integrated light falling 
on that pixel. The signal also integrates a current present in the absence of light, 
the “dark current.” The dark current depends on device temperature. To keep this 
effect constant, and to a minimum, the CCD 1s kept in a vessel at low pressure, and 


is cooled thermoelectrically to —20° C. 


3.2.3 Image Stabilization System 


Compensation is made for small image displacements, caused by seeing and inaccurate 


telescope tracking, by an active tilt mirror in the beam path. The mirror is moved — 


in accordance with signals from a sunspot tracker which senses image motion. The 


active mirror is mounted on three piezoelectric actuators which can tilt the mirror on 


CHAPTER 3. TELESCOPE AND INSTRUMENTATION 44 
lez 


0.8 


0.4 


relative transmission 





0.0 


5894 5895 5896 auth 5898 5899 5900 
MGS 


Figure 3.8: Relative transmission of the SOUP tunable filter vs. tuned wavelength. 
Transmission was measured using a continuum illumination source. 


two axes. The system has a bandwidth of 65 Hz. 

The spot tracker is an array of four photodiodes (a quad-cell). A sunspot image 
centered on the quad-cell will produce equal signals from each diode. If the image 
moves toward one of the four quadrants, however, that quadrant will experience less 
illumination and produce a smaller signal. The differences of signals from the four 
quadrants are used to determine compensatory signals for the active mirror. The 
integrated compensatory signals are used as a feedback to the telescope tracking 


system. 


3.3 Optical Layout 


The optical bench was set up in the Ursies room as shown in the photograph in 
figure 3.9. The schematic lavout of components is shown in figure 3.10. The beam 
from the telescope strikes the active mirror, and is deflected through a field stop 


slightly larger than the active area of the CCD. The heat stop, used to filter out 
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infrared wavelengths, was removed from the hight path because it produced fringing } 
patterns. 

Wavelengths of 5000-7000 A are reflected by the dichroic beam splitter to the. 
tunable filter path. Infrared wavelengths are not reflected. Other wavelengths are) 
reimaged through the Ca Il K-line filter and the 4308 A bandpass filter to the K-lingl 
and “white light” cameras. respectively. The white hght images were used to monitor: 


the seeing conditions and the focus in real-time. The IX\-line camera was not used for 


on a motor driven stage. The magnification of the system is chosen by substituting, 


these observations. 


After reflection from the dichroic beam splitter, the beam is collimated by a lens 


lenses of different focal length. The stage is moved automatically to compensate for} 
the chromatic dispersion of tle telescope. and the whole assembly, moved manually to 
adjust the focus. The beam passes through an aperture stop and is reimaged through 
a set of selectable waveplates and blocking filters, and the tunable filter. For solar 
observations, the beam passes through a sector shutter to the CCD camera. To check 
the position of a spectral line. the flyp mirror was inserted into the beam to deflect it 
to the photometer. 

Part of the beam is diverted. before the polarizer wheel, to the spot tracker. There 
is a neutral density filter and adjustable zoom lens system in this path to optimize 
the intensity and magnification of the image for the spot tracker. 

The motor driven stage. polarizer wheel. blocking filter wheel, and tunable filter 
were controlled by a (PC compatible) COMPAQ computer. The shutter and CCD 
camera were controlled by a VAX Station 3200. The system was free-running, with no 
timebase to control the data acquisition. The timing of the images, therefore, varied 
as shown in figure 3.11. The time mark given to each image came from a crystal 


controlled clock, and is-aceii ate (oa e 
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igure 3.9: Optical bench setup. 
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OPTICAL LAYOUT AT SWEDISH SOLAR OBSERVATORY 
LA PALMA, SUMMER 1990 
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Figure 3.10: 
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Figure 3.11: Timing offset for images vs. time from start of run. The offset is the 
difference between the image time mark and a straight line fit to the time marks. 


3.4 Modulation Transfer Function of the System 


The modulation transfer function (MTF) of an optical system describes the effect of 
aberrations and diffraction on an image formed through the system. If /,(r) is the 
image intensity distribution expected from aberration-free geometrical optics, then 


the real image intensity distribution /,(2) is given by, 
ie) on) NIM ete (a)} S (3.1) 


The attenuation of the spectral power |F{J,}|° is given by the square of the MTF 
functions. 

The components which contributed significantly to the total MTF are the entrance 
pupil, the detector pixel shape, and the earth’s atmosphere. The total MTF is the 
product of the MTF’s of these individual components. The atmospheric MTF will be 


treated in section 5.4.1. 
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3.4.1 Entrance Pupil 


The entrance pupil is defined by the diameter d (48 cm) of the heat shield diaphragm 
in front of the objective. Che MTT of the pupil is given by the autocorrelation of the 


pupil function (a disk). For wavenumbers, 0 < k < kyye, 


We P= 


| to 


lcos7! kK —KVI1 — | (3.20% 


= ( ## pixels k 
where K = ——_—_—__— 


2d \ field of view J byyq 





Bea.e “CCI Camera 


The effect of summing light over the area of a CCD pixel is equivalent to convolving 
the image distribution with the pixel shape, and sampling the new distribution at_ 
points centered on each pixel. Thus. the MTF of the detector is the modulus of the 


Fourier transform of the pixel function (a square). 
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This can be approximated by a cylindrically symmetric function with no more than 
a 2% error. 


MER sinc is (3.4) 


The square of the MTF functions are plotted in figure 3.12 for the entrance pupil and 


CCD detector. 
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Figure 3.12: Squared MTF of telescope and CCD detector 





Chapter 4 
Observations 


The original topic of this dissertation was “sunspot seismology,” that is, the study 
waves propagating into aid out of sunspots. That topic was set aside when waves were 
discovered at substantially higher wavenumber and frequency than had previously 
been observed. The new topic of analysis became the characterization of properties 
of these waves. Many of the observing runs are, therefore, extraneous to the analysis 
described here. They have been included in the list of data in this chapter, but no 


calibration or reduction has been performed on them. 


4.1 Observing Procedure 


The observing runs consist of a repeating sequence of filtergrams with a temporal 
spacing of about 3.5 s between filtergrams. The wavelength, filter width, polariza- 
tion, and relative exposure time for each filtergram in the sequence is defined in the 
sequence file residing on the COMPAQ computer that runs the filter. 

Calibration data was taken before and after observing runs. The order of calibra- 


tion and observations was as follows: 


1. Dark Currents: + measurements of dark currents for each exposure setting. 


Done at the beginning and end of each day. 
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2. Line Centermg: Definition of the nominal line center from which wave- 
length offsets are measured. The photometer detector was used to inte-_ 
grate light over the held of view with the target region centered in the” 
field. | | 

3. Slide Projector Flats: 4 repetitions of flat field measurements using slide 
projector source. 

4. Observing Run: One or more tapes-worth of observations of target region. 


5. Solar Flats: 4 repetitions of flat field measurements using defocused sun- 
light source. 
6, Slide Projector Plitsshenpeat ehcren.3: 
The line centering step is performed by scanning the filter through a range of 


wavelengths (from short to long), and plotting intensity versus wavelength. This step 















was not performed correctly due to a design flaw in the photometer interface with” 
the COMPAQ computer. The computer did not wait for the full signal integration 
time before reading the light intensity. In the resulting line profile scan, the intensity 
points lagged behind the wavelength points. This caused the line to appear shifted 
to longer wavelength. The line was, therefore, “centered” to an erroneously short 
wavelength. This offset of —10 mA for the Na D, line was discovered by making a 


line scan using the CCD camera instead of the photometer (see figure 4.1). 


A? Caltorationw Wetmod 


The purpose of calibration is to produce images with a signal to illumination ratio 
that is independent of both the pixel location and the wavelength setting of the 
tunable filter. The procedure described here is more complicated than the standard 
one used by the Lockheed group. lhis is necessary to provide an accurate calibration 
for dopplergrams. The flat fielding procedure calibrates the whole system, optics and 
detector together. 

The CCD isa very linear device. provided the light intensity stays below saturation 


levels. Response of a pixel can be characterized by a gain and an offset. The offset 
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Figure 4.1: Line scan of Na D, Line using the CCD camera. Signal, integrated over 
the field of view, vs. wavelength offset. 


is determined by the amount of charge that accumulates on the CCD in the absence 
of illumination (the “dark current”). It must be measured for each exposure time. 
Each pixel has a characteristic gain and dark current. The tunable filter is also 
responsible for some variation of the signal strength as a function of wavelength A 
and spatial position z. The net response of the system is to produce filtergrams F;(r) 


at wavelengths A,, given a source intensity distribution S(z, AY 
ERE [S(E)) AZ, A) a + DZ) (4.1) 


where A,(z, A) is the product of filter transmission and CCD gain at filter setting A,, 
and D(z) is the CCD dark current. 
Because the filter bandpass is narrow compared with spectral features in the 


source, we can approximate: 
[SF MED) = ee [Ale ayaa. (4.2) 


Thus 
F(z) = S(#.\,;) B(Z) + D(z) (4.3) 
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whee s Al= / NE igs (4.4) 


In order to determine S, the quantity of interest, we must first measure B; and D. 
The dark current D is easily measured by making filtergrams with no illumination 
source. For § = 0, we have D = F;, 

If we had an illumination source that were both flat in spectrum and flat across 


the field of view, we could measure §; directly. The sources available, however, were: 


1. Defocused Sunlight Ss(A): Flat across field, but not spectrally flat. 


AY 
4 


2. Slide Projector Lamp S;(.t): Flat in spectrum. but not flat across field. 


The technique is to make one filtergram FP at wavelength A, in the continuum } 
(away from the absorption line), using defocused sunlight. And, to make filtergrams | 
FY at \, and F¥ at X, using the slide projector lamp. This gives us the following | 


three measurements: 
Po (®) = Ss(A.) B.(®) + Ds(2) (4.5) 


P(e) = 8, (#) B(®) + Dy (2) (4.6) 
PE(#) = 8, (8) BX) + Dil2) (4.7) | 


Given the dark currents D; and Ds, we can solve for the gain correction G; at A; . 


Gehned sec, tie") ne 


[FY — Dx] 


G(r) = —— 
STE aS | 


(4.8) 
There is some arbitrary overall normalization which depends on Ss(A,). In practice 
the spatial average of G; at the A. filter position was defined as unity, and the other 
gain corrections, normalized accordingly. 
With the gain corrections and dark currents in hand, we can determine arbitrary 


source distributions from filtergrams. 
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4.3 Solar Observations 


The sequence of filtergrams for the sunspot seismology study was named SEISNAFE. 


It was designed to determine the following quantities: 


1. Doppler shift in mid-chromosphere (in sunspot umbrae). 


(H 1 (Balmer) a—line core] 


Le) 


Doppler shift 1n photosphere. 
[Na D,—wings' 


3. Doppler shift in low chromosphere. 
[Na D,;—core] 
4. Doppler shift in photosphere. 


[Fe | (5576)] 


Continuum intensity in green and red. 


[S575A, SS9LA] 


cy 


6. Photospheric magnetic field (line of sight component). 


[Fe I (6302A)] 


The sequence repeated every 17 s. Filter parameters appear in table 4.1. The corre- 
sponding calibration sequences were NAFEFLAT for solar flats, and SEISSP for slide 
projector flats. The observing runs of appreciable length are listed in table 4.2. 

The sequence of filtergrams used for the high-frequency oscillation study was 
named NAFAST. It was designed to make dopplergrams in the low chromosphere 
and photophere, using filtergrams in the core and wings of Na D,, respectively. All 
other measurements were eliminated to increase the repetition rate to once every 
17.7 s. The sequence parameters appear in table 4.3 and the corresponding calibra- 
tion sequences were NAP EF LAT for solar flats, and NAFASTSP for slide projector 
flats. 

The NAFAST sequence was run only once on August 2, 1990 beginning 13:51 
UT. The target region was an active region (AR6181) containing a small, decaying 


sunspot. This is shown on a full-disk magnetogram in figure 4.2. The run lasted 
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Table 4.1: Filter Parameters for SEISNAFE Sequence. 
Tulicere a Polarization 
(0 = linear) | Purpose 


dopplergram 


= “Jt 


continuum 


dopplergram 


wide continuum 


narrow magnetogram 
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Table 4.2: La Palma Observing Runs using SEISNAFE Sequence. 


Date | Time | Duration 
(1990) | (UT) (nun) Seeing FOV | Target Region 


good/v.good | 75” | AR6168 (medium sunspot) 


fair/ var. 150” | AR6168 


























fair/var. 
good/var. 






Variable 






variable 






good/var. 






ValiaDie 













=. good /fair 22: AR6179 (big sunspot) 
46 51 good /fair — overexp'd for umbral osc. 
52 AR6186 (big sunspot on limb) 
quiet region near AR6192 






Table 4.3: Filter Parameters for NAFAST Sequence. 


Offset | Filter Mode | Polarization 
ine (A) (mA) (width) (0 = linear) | Purpose 
95.9 -10 


IE: 
N harrow dopplergram 


harrow ao dopplergram 


5805.9 | 2000] narrow [0 _[ continanim 
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Figure 4.2: Full-disk magnuetogram from the National Solar Observatory at Kitt Peak. 
The field of view containing A\RGIST is outlined. 


L43 min with good/variable seeing that deteriorated in the last 30 min of the run. To 
obtain more information about the active region, the SEISNAFE sequence was run 
before and after the NAF AST runs. .\ continuum image is shown in figure 4.3, a high 


resolution magnetogram, in figure tol. and. an offband H a image. in figure 4.5. 
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Figure 4.3: Continuum intensity (5575A) image of AR6181. The small tick marks 
are oy 
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Figure 4.4: Magnetogram of AR6IS1. Network field is visible around the sunspot. 
The small tick marks are 5”. 
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Figure 4.5: Ha image of AR61I81 (400 mA in the blue wing). The small tick marks 
are 5”. 





Chapter 5 


Data Reduction 


The raw data used for this dissertation reside on a (full) 2.2 Gigabyte Exabyte tape. In 
producing a power spectrum. the volume of data is reduced by a factor of ~ 3 x 10%. 
There is more to the reduction procedure, however, than extracting a manageable 
amount of information from a cumbersome data set. The quality of the data depends 
on careful calibration to account for instrument response, image rotation and drift 
caused by telescope motion. and, distortion caused by the atmosphere. The absolute 
spatial scale of the images must also be determined accurately. This chapter describes 
the details of and cross-checks for the data reduction procedures. 


The reduction steps are listed below in the order they were performed: 


1. Flat Fielding: Calibrate for variation in filter transmission and CCD gain, 
and subtract CCD dark current. 


CCD Flaw Correction: Interpolate over nonfunctional pixels. 


bo 


= 


Correction for Atmospheric Blurring: Normalization of measured MTF. 


ee 


Rigid Rotation: Compensate for image rotation caused by telescope track- 


ing motions during observations. 


Correction for Atmospheric “Shimmer.” (Destretching) 


Doo 


Calculation of Dopplergrams and Line Center Intensity. 


= 


Rigid Alignment: Compensate for image drift caused by telescope tracking 


CLUOrs- 


61 
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S. Fourier Transform of Dopplergrams and Intensities: Described in 


section 2.4. 


9. Angular Scale Determination: Convert units from pixels to arcseconds. 


















| 


; : : : : ; / 
10. Distance Scale Determination: Measurement of foreshortening and direc: | 


tion toward disk center. Described in seétion 6uleIm 


11. Measurement of Oscillation Properties: Described in chapters 6 and 7. 


5.1 Flat Fielding 


The flat fielding scheme is explained in section 4.2. We can use equation 4.9 to 
calibrate filtergrams, given dark currents D and gain corrections G;. The dark current 
is an average of four filtergramms made with a light stop placed in front of the filter. 
The gain correction is constructed using equation 4.8 as follows. 

The three flat field ineasurements F°, F/, and F? are each made from an average 
of four filtergrams. Because the filter stepping motors cid not always find the correct 
positions, it was sometimes necessary to throw out one of the four measurements | 
and use only three to make the average. A filtergram was thrown out if the spatial | 
average of the ceutral 64 x 64 pixels of the image differed by more than 5o from the 
spatial averages of subareas of the other three images, where a is defined by inefficient 


Statistics: 


fmax — !; - 
o, = SX —min (5.1) 
Ny 


The gain correction is the product of (F2> — Ds)7! and (FY — Di) f/f — Dr). | 
The latter is the ratio of the transmittance in continuum to transmittance an Aj. | 
This ratio should be slowly varying over the field of view. The dark currents are | 
not constant with time, however, and the very long exposure times required for slide 
projector flats (2-3 s) caused accumulation of large charges at some small groups of | 
pixels. It was. thus, necessary to “deflaw” F4 and F* (see next section), and smooth | 
the transmittance ratio using a Gaussian point-spread function of 18 pixel fwhm. ; 

The procedure described here produced (flat fielded) images with signal to illu- } 


mination ratios that were constant to within about 3% across the field of view, and, 
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more importantly, to within about 5% across the wavelength range of an absorption 
line. These numbers were measured from filtergrams of defocused sunlight (the flat 


field illumination source) at continuum wavelengths. 


5.2 CCD Flaw Correction 


Of the 2.5 x 10° pixels (in 2 x 2 summing mode), 4513 were not functional. An 
additional 491 pixels had exceptionally high dark currents for exposure times greater 
than about 500 ms. To minimize problems caused by the missing data, values were 
filled in at these pixel locations by linear interpolation from adjacent pixels. This 
procedure is termed “deflawing . 

A good map of nonfunctional pixel locations was available for exposure times 
below 500 ms. For the longer exposure times, a map was constructed using the dark 
current runs. To do this. a cubic trend was fitted to each row of a filtergram. Pixels 
which differed from the trend by more than 7% were flagged as “bad”, and their 


locations, entered in the map. 


5.3 Image Rotation 


The altitude-azimuth axes of mirrors M1 and M2 cause the solar image to rotate 
during the day as the telescope tracks the sun. The rotational direction and rate are 
functions of date and time of day. and can be calculated from geometrical consider- 
ations. Standard software exists to determine the angle of rotation and regrid the 
images by bicubic interpolation. This removes the apparent rotation in a sequence of 


images, but some data is lost from the corners of the images. 


5.4 Atmospheric Distortion 


Turbulence in the earth’s atmosphere will modify wavefronts of light before they reach 
the telescope. Any plane wave passing through the atmosphere will be distorted. 


The resulting wavefront, however, can be approximated by a set of planar patches. 
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The isoplanatic angle describes the angular extent over which image distortions are } 
correlated. For excellent seeing conditions, it can be over 100” (von der Lihe 1988). 
For poor seeing, it can be less than 1”. 

A change in angle of incidence of a planar wavefront produces an image displace- 
ment. When different parts of an image are displaced in different directions, the } 
image appears distorted. The clistortions change with the atmospheric turbulence. 
This can occur over a range of timescales, from | to 1000 ms, with most power in the 


10-100 ms range. 


54.1 “MITE Noriniiahizavion 


The distortion of a wavefront from a given object point differs over the telescope 
aperture. Moreover, local displacements of up to an arcsecond cause the image to 
move around in times of order 10 ms. For exposures lasting tens to hundreds of 
milliseconds, these changing distortions will be superimposed to produce a blurred 
filtergram. The effect is equivalent to convolving the image with some unknown 
point spread function. ‘The Fourier transform of the point spread function is called 
the MTF. It is assumed here that the point spread function is the same everywhere 
in the field of view of a given filtergram. 

The atmospheric MTV can be estimated from the spatial frequency spectrum of 
a filtergram S(f/). A contour plot of spatial frequency spectra is shown in figure 5.1 
for the —40 mA filtergrain set. The sharpest filtergrams in the set are assumed to_ 
be diffraction limited (atmospheric MTF = 1). Since the structure of solar images_ 
should not change appreciably over the observing run, the spatial spectra of diffraction | 
limited images should be a constant Syax(f). Any difference between the spectrum 


of a given filtergram and Sj, must. therefore, be due to the atmospheric MTF. 


Se) = SET ah) : oma (5.2) 


In practice the log-spectra were compared with the mean log-spectrum log S. 


Log-spectra of individual filtergrams were smoothed by convolution with a cosine-bell, 


and the difference between log S and the smoothed spectrum was used as a correctio 
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Log(RMS amplitude) 





O 40 80 120 160 
wavenumber (inverse pixels) 


Figure 5.1: Spatial frequency log-spectra of Na D, —40 mA filtergrams. The solid 
central line is the. mean log-power; the dashed lines are lo contours; the outer dotted 
lines are the maximum and minimum power at a given frequency of any image. 
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factor to the filtergram frequencies. Thus, the "MTF normalization” procedure forced 
the spatial spectrum of each filtergram to look hke the average spectrum. | 

In principle the individual filtergrams could have been normalized to Smax. This 
would have involved very large correction factors for the poorer filtergrams, and would 
have boosted the noise in these filtergrams considerably. Noise in filtergrams produces 


noise in dopplergrams and line center mtensities in a nonlinear manner. So, instead, 


the difference between the mean spectrum and Sa, was retained as a correction 





factor to be applied to the oscillation power. (See chapter 6.) As a final caveat, it 
should be noted that even the best images. assumed above to be diffraction limited, 


actually include some atmospheric blurring. Hence. the correction factor will provide 
| 


y) 


an underestimate for power at the highest wavenumbers. 
The MTF normalization procedure was checked by comparing filtergrams taken_ 
in poor seeing with filtergrams taken at almost the same time in good seeing. After | 


normalization, both filtererams showed the same solar structures. 


5.4.2 Destretching 


They are independent of one another over angular scales greater than the isoplanatic } 


The displacements caused by atinospheric turbulence vary over the field of view. 


angle. Relative displacements between successive images can be measured by local 
correlation tracking (November 1986). This procedure finds the displacement of a 
small subarea of one image with respect to the other from the cross-correlation of the 
two. A displacement vector field is measured using a patchwork of subareas of the 
images. A relative displacement grid of 64 x 64 vectors was made, correlating each 
filtergram with the preceding one. This was done using standard software. 

Local displacements are present not only because of seeing, but also because of 
wave motion and horizontal flows on the sun. We wish to remove the seeing displace- 
ments, but keep the real solar motions. It is possible to distinguish the two on the 
basis of timescale of the motion. Horizontal flows have characteristic timescales of 
hours, oscillations have timescales of minutes, and seeing has a timescale of tens of 
milliseconds. For filtergrams taken seconds apart, seeing motions will appear random, | 


while solar motions will be coherent. 
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The seeing displacement grid was determined from the relative displacement grid 
as follows. For each vector component at each grid location, the relative displacements 
of preceding filtergrams were summed to produce a total net displacement (from the 
first frame). The total displacement was smoothed by convolution with a Gaussian 
of 90 s (5 bin) fwhm, and this smooth trend, subtracted to produce the seeing grid. 

It is important to have high temporal resolution in the data. Correlation tracking 
can lock onto the crests of waves. So, if the impulse response function, used to smooth 
the data, is wider than ~ 120s, the oscillations will be removed from the data along 
with the seeing. The standard Lockheed method for finding the seeing grid is designed 
for horizontal flow stuclies. It uses boxcar averages of 20 or more bins. This would 
be inappropriate for seisniology studies. 

With the seeing displacements measured, the images can be “destretched” to 
their undistorted states. The algorithm to do this treats images as rubber sheets, and 
stretches the sheets (by regridding) based on the local displacement vectors. This part 
of the process also used standard software. The set of filtergrams for each wavelength 
was destretched independently of the other wavelength sets. Destretched filtergrams 
taken at the same time. but different wavelengths, were then visually compared to 


ensure their features lined up. 


5.5 Dopplergrams and Intensities 


Velocity in m/s and line center intensity in arbitrary units were calculated at each 
pixel using equations 2.2 and 2.3. Because of the line centering problem described in 
section 4.1, the line center velocity and intensity were calculated using the —125, —40, 
and +40 mA filtergrams. The actual offsets of these measurements from the average 
line center was about —85.0. and +80 mA. Each dopplergram and intensity image 
was given a time mark averaged from the time marks of the three filtergrams. 

Filter tuning errors aud seeing, which deteriorated toward the end of the run, 
produced some poor quality filtergrams in the data set. This in turned caused some 
nonsensical velocity and intensity determinations. Pixel velocities that were greater 


in magnitude than 10 km/s and intensities that were a factor of 10 greater than 
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average were replaced with the mean velocity/intensity for that frame. 


5.6 Image Drift 

















The rotational regridding assumes the images were rotating about a central point. 
In fact they rotate about the pomt centered in the spot tracker. This inconsistency | 
produces a translational inotion of the solar image. The point (on the sun) centered © 
in the spot tracker, moreover, changes with time. This is caused by variable weighting | 
of parts of the active region as the region evolves and as seeing conditions change. 

From equation 2.27. it is evident that to keep oscillation frequencies from shifting 


by more than 60 plz (0.5 frequency bins) at a wavenumber of 3 Mm7! 


, the image 
drift must be less than 120 m/s. This is equivalent to 1.4 pixel/hr. Image drift is 
measured from the cross-correlation of successive tmages. (It is defined as the distance 
of the maximum of the cross-correlation from the origin.) ; 

The image drift is traced out over time in figure 5.2. It was measured using the 
central 6-4 x 64 pixels of every tenth —125 mA filtergram. As can be seen from the | 
figure, drift can be measured to significantly better than 1 pixel. The measured drift 7 
was smoothed by convolving the time series of each component of the vector with a7 


Gaussian of 6 bin fwhim. The trend was then interpolated to the dopplergram and | 


intensity time marks and removed froin these images by regridding. 


5.7 Angular Scale Determination 


For full disk images of the sun. the image scale can be determined by measuring | 
the radius of the solar disk. For high resolution images, however, the full disk is | 
not visible, and indirect methods must be used. The scale of the La Palma images 
(0.8x lens) was determined by comparison of magnetograms from La Palma and the 
Kitt Peak National Observatory (figures 4.4 and 4.2). The hitt Peak magnetogra 
was made at 15:17 UT on August 2, 1990, and has a scale of 1”/pixel to within 1% 
(Harvey 1991). The La Palma magnetogram was interpolated from magnetograms } 


made at. 13:35 © Panel eee 
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Figure 5.2: Image drift during observing run. Image motion determined by cross- 
correlation is traced out for the observing run. Points represent filtergrams taken 177 
s apart; solid line shows low-pass filter of points; circles represent 18 min time marks 
(t = 0 is at upper right). 


CHAPTER 5. DAZE Diao... 701 










The interpolation method accounted for evolution of the active region by measur- i 
ing the motion of the m. netic field structures using local correlation tracking. The | 
two magnetograms were then stretched by the time-weighted local displacement grid. | 
Finally, the time-weighted signals of the stretched magnetograms were added. 

To compare magnetograms, 43 points on the La Palma magnetogram were flagged | 
manually by cursor on a PINAR computer. Corresponding points were flagged on the } 
Kitt Peak magnetogram. Distances between a point and the centroid of the 43 points 
were then compared for the two magnetograms. The field of view for the La Palma 
images was, thus, determined to be (224.39 + 0.70) Kitt Peak pixels. Assuming the 


Kitt Peak scale measurement. 


La Palma Field of View = 224” + 2” (square) (5.3m 


Chapter 6 


Measurements from Power 


Spectra 


The sets of dopplergrams and intensities were apodized in the space and time di- 
rections with Blackman-IJlarris functions (equation 2.32). Three-dimensional FF'T’s 
were calculated for the full field of view, as well as for the three 79” x 79” subfields 
shown in figure 6.1. The temporal and spatial Nyquist frequencies were 28.3 mHz and 
9.9 Mm’, respectively. [he frequency bin size was 117 Hz, and the wavenumber 
bin size, 0.039 Mm? for the full field, and, 0.11 Mm? for the subfields. 

The squared modulus of each Fourier transform was calculated to produce a 3-d 
power spectrum P(k, vy). The dispersion relations given by equations 1.7 and 1.9 
are symmetric about the 4 = 0 axis. They define “trumpet” shaped surfaces (Hill 


ere (See figure 6.2.) This symmetry allows us 





1988) of approximate form v x 
to integrate the power over the azimuthal coordinate, in order to increase the signal 
to noise ratio. The resultant, collapsed power spectrum P(k,v) will show ridges of 
power at frequencies given by the dispersion relations for the p-modes and f-modes. 

Properties of the oscillations. such as the dispersion relation and mode power, were 
measured from the collapsed power spectra. This was accomplished by fitting one- 
dimensional slices of the spectra. Fits were then compared with theoretical predictions 
to evaluate models of the solar atmosphere. Magnetic versus nonmagnetic regions 


were compared to search for differences in wave properties between the two regions. 


(al 
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Figure 6.1: Magnetogram of ARGISI. The three outlined subfields were used sepa- 
rately to measure wave properties in network field and quiet regions. The small tick 


marks are 9”. 
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Figure 6.2: Surfaces of power in (4, v) space. 
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6.1 Collapsing 3-d Power Spectra 




















The azimuthal symmetry of equations 1.7 and 1.9 is realized only in solar coordinates, | 
comoving with the gas. Iu practice, one or both of these conditions may not be met. 
The frequency shift described by equation 2.27 occurs, with a cosine dependence on 
azimuth @, if there is a relative (horizontal) motion between the observer’s coordinates | 
and the gas. The rigid alignment procedure described in section 5.6 keeps the observa- 
tional coordinates fixed with respect to the sunspot. It has been shown, however, that 
there is relative motion between magnetic flux and nonmagnetized gas (Snodgrass & © 
Ulrich 1990). Hill (1988) fines flows of ~ 100 m/s below the photosphere. 

Visual inspection of the power spectra did not reveal evidence of flows. Thus, . 
if present, flows produce frequency shifts of less than one frequency bin (117 «Hz). . 
If plane waves propagated isotropically. a sinusoidal frequency shift with @ would | 
serve to widen the ridges in the collapsed power spectra. In the measured spectra, | 
however. power depends mildly on angle 8. Thus, a horizontal flow would produce 
slight asymmetry in addition to broadening a ridge. 

The second condition affecting the symmetry of the trumpets is geometrical. Fo 
observations away from clisk-center. the projection of solar coordinates onto obser 
vational coordinates causes foreshortening. Waves propagating toward or away fro 
disk-center appear with shorter wavelength (higher wavenumber) than waves prop- 
agating in the perpendicular direction. This foreshortening distorts the cylindrical 
symmetry of the trumpets to give them an elliptical rather than circular cross section. 


Procedures to measure aud compensate for this effect are described below. 


6.1.1 Distance Scale Determination 


The distance scale of the images or, equivalently, the wavenumber scale of the power | 
spectra follows from trigonometry. once the magnitude and direction of foreshortening | 


have been determined. This measurement was made by fitting an ellipse to the cross 


7 


is shown in figure 6.3. The slice has been remapped to cylindrical polar coordinates. | 
' 
| 
| 


section of a trumpet in a 3-d power spectrum. 


A constant frequency slice. in the 5-minute band, of the full-field Doppler spectru 


==] 
Vit 
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Pee 6.3: Constanteiveqtcweresitce of thes3-depower spectrum. A shce at 3.99 mHz 
has been remapped to cylindrical polar coordinates. Log( power) 1s shown as a function 
of wavenumber and azimuth. 


So, the cross section of a trumpet appears at roughly constant wavenumber for all 
values of 8. These horizontal bands of power are the “rings” described by Hill (1988). 


In polar coordinates. an ellipse centered at the origin has the form, 


|— 


= eel cos (0 0o)]7 (Oa 


where 49 1s half the minor aNnis (the true wavenumber); 99 is the angle of the major 
axis (direction toward disk-center): and sis the ratio of minor to major axes (the 
fractional foreshortening). A close inspection of figure 6.3 shows that each of the 
bands exhibits this oscillatory behavior of k as a function of @. 

The values of fy. 9. and yp were fitted to the measured wavenumber /(@) for the 
Pp; ring in the frequency slice vy = 3.99 mHz. The function 4(@) was defined as the 


wavenumber at peak power in the constant frequency slice. It was determined by 
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quadratic interpolation. [he measured fractional foreshortening was, 


it = 09au 














This corresponds to a distance from disk-center of 17°, which agrees with an ap prox: 
imate measurement made from the full-disk magnetogram in figure 4.2. 

The assumption of a planar atmosphere has been made here, in treating the fore 
shortening as constant across the field of view. If the foreshortening varies appreciably, | 
modes of different wavenumber will be mixed in the 3-d power spectrum. This will 
smear the trumpet surfaces in (he. v) space in a manner which cannot be corrected. : 
For large fields of view. the images must be mapped to “flattened” coordinates before 
transforming to Fourier space. The aplanarity for the data analyzed here was not | 
great enough: to warrant such treatment. } 

The 22-4” fleld of view represents 13° along a great circle on the Sun. Since the q 
average distance from disk-center was about 17°, plane waves propagating along the 
direction of foreshorteniny will be mixed over a wavenumber range of +3.5%. While 
this may seem large, only the waves propagating exactly along a diameter of the disk 
experience mixing to this extent. Waves propagating in the perpendicular direction 
are not mixed at all. The net widening of the ridges is, thus, small. This can be seen 


in figure 6.3. The width of the ridges appears independent of 8. The width, therefore, 


is not an artifact of the coordinate system. 


6.1.2 Azimuthal Projection 


The power was integrated around ellipses in the (A;, fy) plane. To do this, ko was 
computed from equation 6.1 for each ordered pair (A,, ky) in the 3-d spectrum, and 
power from that bin was added to the appropriate & bin in the 2-d spectrum. Bins 
in the 2-d spectrum were normalized by (multiplying by) the ratio of the area of an 
elliptical band to the number of bins accumulated from the 3-d spectrum. The area 


A of an elliptical band was approximated by, 


fou-Ve+P? (6.: 


Iw | 3 
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Figure 6.4: Scaled velocity-power spectrum. A portion of the collapsed power spec- 
trum for the full field of view. with variable gray-scaling in the & and v directions. 
Background has been subtracted. The noise visible at the top of the spectrum is an 
artifact of the scaling. From Fernandes ¢f al. (1992). 


where a and ) are the major and minor axes of the ellipse; and w is the width of a 
wavenumber bin in the 2-d spectrum. The two-dimensional spectrum P(k,v) for the 


full-field dopplergrams is shown in fi@ure 6-4. 
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6.2 Parametrization and Least-Squares Fitting 


6.2.1” Fitting Procedure 


The frequencies, widths. and heights of the ridges were measured as functions of 


vidually fitted by least-squares method. [Examples of slices at constant wavenumber 


| 


appear in figure 6.5. A single isolated mode should have a Lorentzian profile. 


Palen 


Jes (6.3) 


Where / is the mode power: vp 1s the mode frequency; and [ ts the fwhm. The high 
frequency mock modes described by Kumar eft al. (1990) should have a sinusoidal 
profile. 


sin’(watgy ) 


Pe (6.4) 


wavenumber. To do this. constant wavenumber slices of the power spectra were indi- | 
w? — w,* sin*(wa/c) 

where 2g 1s the source depth: «, 1s the acoustic cutoff; and a is the height of the upper 9 
reflecting surface. Solar processes such as horizontal flow fields, and measurement 


effects such as seeing and apodization tend to smear the shape to produce a Gaussian 





profile. 
i  — fy tase (6. 
Tey 2 ' 
where Py is the mode power: vy is the mode frequency; and 0 = fwhm/2.36. 
The ridges he on top of a background which decreases monotonically with wave- 
number. In the photosphere, the granulation is caused by overshoot of convective cells 


into the overlving convectively stable lavers. The power, in the photosphere, below 






: 


the f-modes is caused by this convective overshoot. Extrapolation of calculations of 





Nesis & Mattig (1959) to the low chromosphere suggests a similar explanation for the 
background power there. 

Because of the uncertainty of the physical description for the background, ther 
is no proper way to parametrize it. Nevertheless, the background must be accounte 
for to avoid contaminating the ridge fits. [t can either be included as an extra term 1 
the fits, or fit independently and subtracted from the power spectrum before fittin 


the ridges. 
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Figure 6.5: Constant wavenumber slices of the 3-d power spectrum. Data are plotted 
as points with error bars. Background and individual Gaussian modes are plotted 
as dashed lines; and the sum of background and modes is plotted as a solid line. 
Top left: & = 1.16 Mm’; Bottom left: & = 5.25 Mm7’; Right (top & bottom): 
>= 2.34 Mm. 
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The latter method was chosen because the background could be well characterized 
using the range of frequencies on either side of the ridges. The shape of the background 


P, suggested a power law parametrization. 
leur) = eve (6.6) 


where A, and By, are parameters fitted independently for each wavenumber slice. 
The frequency ranges used for the background fits include: (0.12 < v/mHz < 
1.3) and (12.0 < v/mHz < 15.0). 

The ridges were fitted after subtracting the background. As can be seen from 
figure 6.5, the ridges overlap. and so must be fitted simultaneously. With three 
parameters per Gaussian. and up to eight visible ridges, the fits involved up to 24 pa- 
rameters at a time. This large parameter space combined with noisy data necessitated 
the development of a very robust fitting algorithm. 

The fitting program is described in the next section. It provides a nonlinear, 
weighted, least-squares fit to ungridded data. The weights used in the least-squares 
fit were the measurement errors. Hence, the fitting procedure minimized y*. The 
error on the power was estimated from equation 2.38. This in turn assumed the 


velocity and intensity errors derived in section 2.1. 


6.2.2 Levenberg-Marquardt Method for Nonlinear Regres- 


sion 


The basic algorithm follows the Levenberg-Marquardt procedure (Press et al. 1989). 
The method attempts to minimize the \~ in parameter space (@), given N measured 


Valles iy O, al poltser 


eres on ar 2 
ela = y (4 Jes 2) (6.7) 


where y(v:d) is the fitted function. 


The minimum of the surface \7(@) is searched for in an iterative manner. Succes- 


sive guesses for the parameter vector @ are made based on the gradient and curvature 





| 
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matrix for the \? surface. When the distance to the minimum, as estimated from 
the curvature, is large, the gradient is weighted more heavily in the guess. When the 
minimum is near. the curvature carries more weight. The algorithm assumes it has 
found the minimum when it can no longer reduce \? by more than 0.1. The covariance 
matrix for the fit is simply the inverse of the curvature matrix at the minimum. 

To prevent runaivays. the parameters are tested after each iteration of the algo- 
rithm in which \* is reduced. Parameters that are not within set limits are fixed at 
appropriate values, and no longer varied in the fit. The remaining parameters are 
reset to their initial valucs @. and the procedure starts again from the beginning. 

For Gaussian ridges, the limits for the height parameter are € (small, positive) 
and Imax (very large). Values falling outside the limits are fixed at the limit. The 
limits for the width parameter are one-half and twice the expected width. Again, 
values falling outside the limits are fixed at the limit. The limits for the frequency 
parameter fall halfway between the expected positions for ridges (a). For the bottom 
(top) ridge in a fit. the lower (upper) limit is defined such that the two limits are 
equidistant from the expected ridge position. Values falling outside the limits are 
fixed at the initial guess (o. 

The reason for this treatment of parameter limits is to prevent noisy ridges from 
disturbing the fits of adjacent ridges. Ridges which are too noisy to be properly fitted 
are treated as background. 

For the first spectrum to be fitted, initial guesses at a number of wavenumbers were 
made by eye, using cross hairs on a workstation display. Guesses at other wavenum- 
bers were interpolated or extrapolated from the nearest four fitted wavenumber slices. 
For the remaining power spectra. initial guesses were interpolated from the already 
fitted spectra. 

In addition to the background and ridge fits, the least-squares fitting program was 
used to interpolate and extrapolate mode frequencies to make the initial guesses, and 


to determine the parameters of equation 6.1 in section 6.1.1. 
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Figure 6.6: Power law exponent for background fit. Exponent vs. wavenumber. 


6.3 Results of Fits 


6.3.1 Background 


The exponent to the power law fit 1s plotted in figure 6.6. The background velocity 
power, integrated over frequency, is plotted in figure 6.7. The spectrum falls off much 
more slowly with wavenumber than the spectrum of convection in the photosphere re- 
ported by Chou et al. (1991). The root-mean-square (RMS) velocity is about 820 m/s. 
Nesis & Mattig (1989) have measured RMS velocities, due to “mechanical overshoot”, 
that increase with height in the upper photosphere. Their model predicts the height 
dependence, for line-of-sight velocities, shown in figure 6.8. A smooth extrapolation 
of this curve to the height of formation of the Na D, line (500-700 km) yields an 


RMS velocity of 700-1100 m/s, in agreement with our measurements. 
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Figure 6.7: Background velocity power integrated over frequency. Background was 
fitted from the “quiet” subfield spectrum. From Fernandes et al. (1992). 


6.3.2 Mode Frequencies 


Frequencies of ridges froin the full-fleld velocity and intensity spectra, and the quiet 
subfield velocity spectrum are tabulated with their respective statistical errors in 
appendix A. The velocity points are plotted in figure 6.9. All tabulated points 
correspond to fits in which the offset. width, and height of the Gaussian converged 
to physically reasonable values. and the integrated power under the Gaussian was at 
least 4 standard deviations above the background. 

All six fitted p-mode ridges extend up to at least 7.0 mHz. The ps and pe ridges 
extend up to a frequency of 9.5 mHz. There is no deviation from a smooth extrap- 
olation of the ridges at lower frequencies. Libbrecht, Woodard & Kaufman (1990) 
have made measurements that overlap with ours at frequencies up to 5.3 mHz. As 
can be seen from figure 6.10. the fitted frequencies are consistently higher than those 
of Libbrecht et al. by an average of (35 + 13) wHz. This difference, however, is only 


a quarter of the frequency bin width, and is indicative of the systematic error in the 
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Figure 6.8: Calculated Doppler velocity vs. height for the Fe I (5123.75 A) line. From 
Nesis & Mattig (1989). 


measurement. 


The f-ridge is visible at wavenumbers up to 5.57 Mm7’. 


There is, however, a 
systematic discrepancy between frequencies fitted from the velocity and intensity 
spectra as can be seen from figure 6.11. Frequencies measured from the velocity 
spectra are consistently higher than those from the intensity spectra by up to 50 wHz 
at wavenumbers below 2.5 Mm7!. This difference increases to ~ 100 Hz at 4 Mm™". 
The discrepancy may in part be due to a difference between background shapes in 
the velocity and intensity spectra. 

Libbrecht et al. (1990) report a Gaussian background, centered near 3 mHz, in 
their velocity spectrum. Inclusion of such a background in our velocity fit would bring 
the f-ridge frequencies into better agreement with the intensity data at wavenumbers 
below 1.1 Mm7', but would push them farther apart at higher wavenumbers. Fitting 
such a background function, moreover, would require a longer time series with corre- 
spondingly better signal-to-noise ratio. It is noteworthy that the f-mode frequencies 


measured from the intensity spectra differ with those reported by Libbrecht et al. by 


CHAPTER 6. MEASUREMENTS FROM POWER SPECTRA 85 


y (mHz) 





—0.4 0.0 0.4 0.8 
Log(k / Mm7') 


Figure 6.9: Measured mode frequencies (from velocity spectra): v vs. logig’. Fre- 
quencies of the fand p,—p, ridges are plotted with 1 o error bars. For k < 3.82 Mm", 
frequencies are taken from the full-field spectrum. At higher k, values from the quiet 
subfield are plotted. 
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Figure 6.10: Comparison of measured p-mode frequencies with Libbrecht et al. (1990). 


Frequency difference vs. & is plotted for p,, p2, p3 (from top down). 
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Figure 6.11: Comparison of f-mode frequencies measured from velocity and intensity. 
Frequency difference vs. -.. 


under 30 «Hz, as shown in figure 6.12. 

The weighting of magnetic regions is different for intensity and velocity spectra. 
However, as shown in section 6.3.5, the discrepancy is not caused by a difference in 
mode frequencies in the magnetic and nonmagnetic regions within the field of view. 
Regardless of which set of f-mode frequencies we use, the measured frequencies lie 
below the theoretical frequencies given by equation 1.9 at wavenumbers greater than 
1 Mm7’. The difference is 100-200 ;.Hz at 4 Mm7'. The fractional difference between 


measured and theoretical frequencies is plotted in figure 6.13. 


6.3.3 Mode Widths 


Ridge widths are plotted as error bars to the frequencies in figure 6.14. The ridge 
widths at low frequency are due to the apodizing function. At high frequency, the 
ridges are broader than expected for mock modes. The measured widths probably do 


not represent the true widths. The width is the most sensitive parameter to noise in 
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Figure 6.12: Comparison of f-mode frequencies from full-field intensity spectrum with 


Libbrecht et al. (1990). 


the fitted power spectrum. 


6.3.4 Mode Power 


Mode power is defined as the integral under the fitted Gaussian. 


Y2r 


P= 
2.36 





hw (6.8) 


where A is the height of the ridge; and w, the fwhm. The error on the power ts given 
by, 

op = Py (28) + (2) 2m , 
where np," is the covariance of the height and width. The mode power is plotted 
as a function of frequency in figure 6.15. The p-mode power has the same frequency 
dependence for all p-ridges. There is very little explicit wavenumber dependence. 


Power increases as v'* to a maximum at 3.4 mHz, at which point it turns over and 
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Figure 6.13: Comparison of measured f-mode frequencies with theory. Frequencies 
from the full field intensity spectrum are compared with the relation w*? = gk (g = 
eee m/s’). 
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Figure 6.14: Frequency and ridge width vs. &. Top plot shows f-modes and even 


order p-modes. Bottom plot shows odd ordered p-modes. Error bars represent the 


fwhm from the fit. 
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Figure 6.15: Velocity power vs. mode frequency. From Fernandes et al. (1992). 


falls off again as y~*° 


. The f-modes have less power at a given frequency than the 
p-modes. As surface gravity waves with very little compression they are physically 
different from the p-modes and are not expected to show the same power spectrum. 

An MTF correction was made to the power plotted in figures 6.7 and 6.15 by com- 
bining the calculated MTF of the telescope and CCD with the estimated atmospheric 


MTF. See sections 3.4 and 5.4.1. 


6.3.5 Comparison of Magnetic and Nonmagnetic Regions 


Mode frequencies were compared for the quiet, canopy, and sunspot subareas shown 
in figure 6.1. The frequency differences are plotted in figures 6.16 and 6.17. The 


frequencies of ridges in.each subarea agree with those of the other subareas to within 
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Figure 6.16: Mode frequency differences for canopy and quiet subfields. Top left, 
f-mode; top right, pj; bottom left, p2; bottom right, ps. 


the scatter of the data (+40 ,Hz). 

The experimental scatter is much larger than the error bars derived from the fit. 
This may simply indicate an underestimation of the velocity error, or it may be a 
result of solar noise. Modes of finite lifetime have an intrinsic frequency uncertainty. 
Determination of the frequency of the cavity resonance vo requires observations for 
several mode lifetimes. The 2.4 hr observing run probably constitutes less than one 
mode lifetime over much of the wavenumber range. Hence, the measured frequencies 
will be distributed about vo in a manner which depends on the mode lifetime, and 


shape of the wavepackets. 
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Figure 6.17: Mode frequency differences for sunspot and quiet subfields. Top left, 
f-mode; top right, p;; bottom left, p2; bottom mght, ps. 





Chapter 7 
Spatial Distribution of Power 


The spatial distribution of power was measured for groups of modes, as well as for the 
sub-fundamental region of the power spectrum. To do this, the set of dopplergrams 
V(z,t) was masked using W(.v) from equation 2.33, and a 3-d FFT, calculated for the 
full field of view. A spectral mask TWEE. v) was then applied to the Fourier transform 
to remove all but a small range of wavenumbers of a given ridge. An inverse Fourier 
transform of this masked 3-d spectrum yielded velocity Wy,;(%,t) for a specific set of 
modes. 


Wasa. t) = Fo) Mh) PC (az) (ee. t)}} (7.1) 


The mean power was calculated at each pixel to make a velocity power map Py;(2). 


* eet: = 
Py (2X) = =f pel eeclliie dt (ie) 
| "Vnygq ~ 
7 - | FSM Ry) FW (@)V(B.O}}] dv (7.3) 


where F5' denotes the spatial, inverse Fourier transform. 


7.1 Mode Masks 


Masks were constructed for the f through p3-modes in 6 wavenumber ranges. In addi- 
tion, two sub-fundamental and a high frequency set of masks were made for the same 


wavenumber ranges. The mask limits are shown in figure 7.1. Three-dimensional 


9.4 
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masks were made from tliese 2-d limits by rotating the loci in the azimuthal direc-7 
tion. To account for the foreshortening. the rotation followed ellipses described by 


equation 6.1. 


7.2 Results 


Two sample power maps are compared with a magnetogram and an H @ image in 
figure 7.2. Only the central disks of the frames are visible due to the apodizing 
window. The sub,-f and f-inode maps were both made for the wavenumber range} 
1.95 < k/Mim7™! < 2.34. This is equivalent to a wavelength of about 4”, or 3¥imm 
The mottled appearance of the power maps is an artifact of random interference of” 
waves in a finite time series. ‘The sub,-f map shows little variation other than an] 
overall trend of higher power in the upper right hand corner. The f-mode map shows 
a depression mm power in the active region. | 

A selection of power maps is shown in figure 7.3. An estimate of the RMS velocity 
i each map can be obtained from figure 7.4. The power in high wavenumber modes 
is depressed in the active region. At low wavenumbers, the active region has little 
effect on the power. The sub,-f and sub2-f maps show no correlation of power with 
magnetic field at any wavenumber. The high frequency maps, though not shown in 
figure 7.3, looked very siinilar to p-mode maps of the same wavenumber range at 
lower frequency. There is little difference between the f, p;, pz, and p3 maps for a 
given wavenumber range. 

A horizontal slice through the />mode map from figure 7.2 was made above the 
sunspot. The slice was low pass filtered using a Gaussian impulse response function © 
(16 bin fwhm). See figure 7.5. This low pass filtered slice is compared with the 


magnetogram signal in figure 7.6. The power depression is significantly broader than 


ry e ry » | 
the spatial extent of the photospheric magnetic field. It is narrower than the canopy 


field at ~ 1200 km. as can be seen from the H @ image in figure 7.2. The filaments 
visible in H a are streams of gas traveling along magnetic field lines. The filaments, 
thus, trace out the field lines which form a canopy above the photosphere and lower 


chromosphere. 
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Figure 7.1: Masks for power distribution maps. Measured mode frequencies are 
shown as dots. Mask limits are shown as solid lines. The high frequency (HF') masks 
covered the range 6.0 < vy/mHz < 10.0; the p-mode mask limits were set halfway 
between modes; the f-mode limits were set equidistant from the f-mode; the sub2-f 
band extends from the f-mode down to v = 0; the sub,-f band is the same width as 
the f-mode band. It is a subset of the sub2-f band. 


CHAPTER 7. SPATIAL DISTRIBIGFION OF Pow? 97 


2 ddd geht ike ow ote me ili mh 


howe og ww a 


4 
J 


* 
Oe ee ee Re Te Oe ee ee ee 





ligure 7.2: Comparison of power maps with maenetic fields. Top left: log power 


. ° ae rs | ; . ° . 
for sub,-f (1.95 < F/Nimn 2.51). Top mght: absolute value of magnetogram (log 
scale). Bottom left: Ha — 100 mA image. Bottom right: log power tor f-mode 


- : / —1 porte Nae : : - 
(1.95 < A/Min < 2.3-1). The small tick marks are 5”. 
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Figure 7.3: Power distribution for sclected masks. Log power is displayed for sub.-f, 
f. Pr. pg. and pg maps with wavenumber ranges: 
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Figure 7.4: RMS velocity vs. &. Curves fitted from quiet subfield velocity spectrum. 
Curve labeled sub-fis taken from the background fit. 


It is surprising that the sub-fundamental power is not correlated with magnetic 
field. The RMS convective power in the photosphere is a factor of 2 lower in plage 
(network field) than in nearby quiet regions (Title et al. 1992). Moreover, the sub2-f 
band hes just below the /-mode band, but shows completely different behavior. At 
the same wavenumber, and almost the same frequency, oscillation amplitudes are 
reduced in active regions. whereas non-oscillatory motion is not. 

It is impossible to distinguish between instrumental backgrounds and solar back- 
grounds from a single data set. It is difficult, however, to come up with an obser- 
vational effect that would produce the measured background. Filter transmission 
variation should affect only the very low wavenumber part of the spectrum. Seeing 


variations should be independent of frequency; whereas, the observed background 


decreases with frequency. 
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Figure 7.5: Slice of f-mode map. Horizontal slice from the (1.95, 2.34] Mm7' map. 
Also shown is a low pass filter of the slice. 
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Figure 7.6: Comparison of f-mode slice with magnetic field. Dashed line shows low 


pass filter of slice (arbitrary scale). Solid line shows magnetogram signal (arbitrary 
scale). 


Chapter 8 


Conclusions and Discussion 


8.1 Surfaces and Cavities 


The existence of a chromospheric cavity 1s expected to produce a resonance consist- 
ing of a horizontal band of power and “avoided crossings” at about 5.5 mHz. (See 
figure 1.4.) This cavity is required by the model of Balmforth & Gough (1990) to 
explain the existence of high frequency modes. The p-mode and f-mode ridges shown 
in figures 6.4 and 6.9 are clearly visible at 5.5 mHz, but there is not a horizontal band 
of power nor are there “avoided crossings. The clear absence of this resonance in our 
data shows that there is no chromospheric cavity. Either the reflection coefficient at 
the base of the corona is too low or the surface is too irregular to form a good cavity. 
This result supports the alternate hypothesis of Kumar et al. (1990), that the 
high frequency ridges are a purely geometrical artifact. The result is in agreement 
with the conclusion of Woodard & Libbrecht (1991): and it contradicts the claim of 
Dame et al. (198-) that there is a horizontal ridge of power at 5.5 mHz. Our power 
spectrum has a significantly better signal to noise ratio than that of Dame et al. 
The f-modes are sensitive to conditions in the chromosphere and transition re- 
gion. The fact that f-mode frequencies are lowered at high wavenumber may indicate 
interaction with chromospheric magnetic fields. energy leakage into the corona, or 


horizontal inhomogeneitics in the chromosphere-corona boundary. The last of these 
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three possibilities is interesting from the standpoint of acoustic cavities. Perhaps ir 
regularities in the interface between the chromosphere and corona are responsible for 
both lowering the f-mode frequencies, and for reducing coherent reflection of acoustic — 


Waves. 


8.2 Waves in Magnetic Fields 


There is no frequency shift for f and p-modes in magnetic network fields of 100 — 
—1000 gauss to within the scatter of the data (~ 50 «Hz). This limit is not stringent | 


enough to determine whether the frequency shifts measured by Woodard et al. (1991) } 
























result from small shifts over the whole disk. or larger local shifts, in active regions, | 
which contaminate the full-disk measurement. The modes measured by Woodard et. 
al. have widths of order | jfIz. Hence. to contaminate such a measurement, mode 
frequencies in active regions would have to be distributed within a few pHz of the 
quiet sun values. A measurement of such precision would require observations of at. 
least 20 hr duration, assuming the waves are coherent for the observing period. At 
La Palma, 20 hours of observations can be made in two davs. Hence, the waves would | 
add coherently if the mode lifetime is > 34 hr. [f the mode lifetime is less than this, 
the signal to noise ratio increases only as the square root of time, so a much longer 
observing time would be necessary. 

Power maps for oscillations in the low chromosphere show that, at high wavenum- 
ber, wave amplitudes are suppressed 15 Mm outside of the photospheric magnetic | 
field. (See figure 7.6.) In the photosphere. the power is suppressed within one (hort- 
zontal) wavelength from the field (Tarbell ¢¢ af. 1988). The vertical distance between 
the photospheric and chromospheric observations is under 0.56 Mm, and the (true) } 
wavelength of the oscillations is ~ 2 Mim. Diffraction cannot account for the spread- 
ing of the region of the power depression because the waves have traveled less than 
one wavelength in the interval over which the spreading occurs. The power suppres- 
sion must, therefore, be produced by interaction of the oscillations with the field zn 
the layer in where the waves are observed, 


The area containing magnetic flux was not measured directly at the height of the 
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observations. Its extent can be inferred, however, from the magnetogram made in 
the mid-photosphere, and the H a image in the mid-chromosphere. The field lines 
spread out rapidly above the photosphere. Hence, the extent of the field at the 
observed height should be intermediate between the limits of the mid-photospheric 
and mid-chromospheric fields. This is consistent with the above conclusion. 

The f-modes propagate without compression, so they produce no pressure fluctu- 
ations in the region where oscillations are suppressed. The mechanism suppressing 


the wave motion must, therefore, act through the velocity fluctuations. 


8.3 Directions for Future Work 


The frequency measurements presented here show the existence of high frequency and 
high wavenumber oscillations. They show general characteristics of the dispersion 
relations such as the smooth continuation above the acoustic cutoff. They are not, 
however, high precision measurements in comparison with existing low wavenumber 
observations. Precision of the measurements can be improved with longer observing 
times. This would improve the signal to noise ratio, and enable better fitting of the 
background and ridges. With better frequency resolution, it would become necessary 
to correct the fitted frequencies for their weighting with mode power (Ixorzennik 
P91). 

Important characteristics of the solar oscillations can be inferred from their phases 
(Deubner & Fleck 1989, 1990: Fleck & Deubner 1989). Trapped and nontrapped 
oscillations can be distinguished by the phase difference between velocity and intensity 
measurements. Propagating and evanescent waves can be distinguished by the phase 
lag between velocities at different heights in the atmosphere. An attempt was made 
to measure velocity-intensity phase differences. The attempt failed, possibly because 
of the large height range included in the measurement (see section 1.4), and possibly 
because of the large number of modes included in a single frequency bin (Braun et 
al, 1992a). 

The interaction of magnetic fields with oscillations can be better quantified by 


measuring both wave amplitudes aud magnetic fields at a range of heights in the 
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atmosphere. Improvements that have been made to the tunable filter system since 
the 1990 observing seasou will allow Doppler measurements restricted to either the 
core or wings of the Na D, line. In addition, magnetograms can be made in the core 
or wings of the line by finding the difference between dopplergrams made with left 
and right circularly polarized light. 

Braun et al. (1992b) report halos of enhanced high frequency acoustic power 
around active regions. A search was made for such halos in power maps made with 
the HF inasks shown in figure 7.1. No enhancement of power was found. This does 
not contradict the result of Braun ef af. as they measured intensity fluctuations in 
the photosphere. in a lower wavenumber regime. It would be interesting to repeat 
the analysis with continuum brightness measurements for a field of view of at least 
6 arc minutes. 

The cause of the sub-funcdamental background is, currently, unknown. Experi 
ments should be done to discriminate between observational and solar backgrounds. | 
For example, a set of filtergrams taken at fixed wavelength with a high cadence 
could be used to make artificial dopplergrams, and check the effect of seeing on back- 
ground. Another set of fillergrams taken at the nominal, observing wavelengths could 
be made, using integrated sunlight, to check the effect of filter-transmission variation 
on background. 

An attempt was made to find the oscillation sources described by Brown (1991). 
temporally high-pass filtered movie of velocities was scanned visually for evidence of | 
spatially and temporally localized sources of power. No good candidates were found. 
A better method might be to construct the two-dimensional acoustic Pointing vector 
described by Braun ef al. (1992b). Large positive divergences of this vector field for 


high-pass filtered data could point out oscillation sources. 


Appendix A 


Frequency ‘Tables 


Table A.1: Mode frequency v +0, (Hz) as a function of wavenumber k (Mm7') — 
measured from dopplergrams. For & < 3.82 Mm7', frequencies are taken from the 
full-field spectrum. At higher &. values from the quiet-subfield are tabulated. 
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Table A.2: f-mode frequency v + 0, (Hz) as a function of wavenumber k (Mm7') 
— measured from intensitles. 





APPENDIX A. FREQUENCY ing. 111 








32995 326032 3.669 3.707 3.745 3.783 
Le 4874 15 4910 17 4914 16 4948 16 4896 20 4961 20 
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